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I.  Molecular  Genetic  Studies  of  Bone  Mechanical  Strain 

1 .  Molecular  Genetic  Studies  of  Bone  Mechanical  Strain  -  in  vivo  studies 
Introduction 

Physical  exercise  is  believed  to  increase  bone  formation  and  improve  bone  density  in 
osteoporotic  patients  and  other  patients  with  metabolic  disease  as  well  as  in  novel  subjects. 
However,  the  mechanisms  by  which  mechanical  stress  converts  mechanical  signaling  and 
activates  a  set  of  stress-responsible  genes  leading  to  bone  remodeling  are  not  fully  understood. 
Our  previous  studies  have  shown  that  two  inbred  mouse  strains,  C57BL/6J  (B6)  and  C3H/HeJ 
(C3H),  respond  differently  to  mechanical  loading  in  vivo.  B6  mice  showed  a  significant  increase 
in  bone  density  as  compared  to  C3H  mice  when  subjected  to  similar  loading  regimens.  These 
results  suggest  that  mechanical  stimulation  may  induce  distinct  types  of  gene  expression  in  B6 
mice,  which  contributes  to  the  subsequent  bone  formation.  Therefore,  it  is  of  particular  interest  to 
us  to  localize  the  chromosomal  regions  responsible  to  mechanical  stress  and  identify  osteogenic 
candidate  genes  by  applying  powerful  genetic  tools,  namely  quantitative  trait  loci  (QTL). 

Technical  Objectives 

Technical  Objectives  1  and  2  for  the  past  twelve  months  was  revised  on  04/15/03.  The 
revised  specific  objectives  for  Technical  Objective  1  are  provided  below.  In  the  original  proposal, 
there  were  two  Technical  Objectives  for  this  portion  of  our  work.  The  second  Technical 
Objective  has  now  been  incorporated  into  the  first  Technical  Objective.  Therefore,  on  this 
progress  report  we  will  report  on  Technical  Objective  1,  and  its  five  specific  objectives.  With  the 
aforementioned  revision,  there  is  no  longer  a  Technical  Objective  2;  Technical  Objective  1  will 
then  be  followed  by  Technical  Objective  3,  which  is  a  progress  report  on  our  in  vitro  work,  to 
evaluate  the  physical  actions  of  mechanical  strain  on  bone  cells  in  vitro. 

Our  specific  objectives  for  the  revised  Technical  Objective  are: 

1)  To  establish  the  optimal  method  for  inducing  mechanical  strain  on  the  tibia. 

2)  To  determine  the  optimal  method  for  quantitating  the  bone  formation  response  to 
mechanical  strain  on  the  tibia  in  vivo. 

3)  To  determine  the  number  of  days  required  of  4-point  bending  induced  mechanical 
strain  on  the  tibia  in  order  to  elicit  an  optimal  response  for  quantitative 
measurements. 

4)  To  select  the  optimal  mouse  inbred  strain  pair  (i.e.  B6  and  C3H  or  some  other  strain 
pair)  to  perform  the  QTL  mapping  study. 

5)  To  determine  the  appropriate  age  of  animals  to  apply  the  4-point  bending  in 
order  to  obtain  valid  quantitative  measurements  of  the  bone  formation 
response. 

Our  goal  is  to  evaluate  the  bone  formation  in  1 0  different  mouse  strains  using  mechanical 
loading  in  order  to  identify  two  mouse  strains  which  show  extreme  differences  in  the  bone 
formation.  These  selected  strains  will  then  be  backcrossed  for  two  more  generations  (F2)  and  the 
F2  mice  will  be  used  to  identify  candidate  genes  responsible  for  bone  formation  by  applying  a 
powerful  genetic  tool,  quantitative  trait  loci  (QTL),  as  described  previously.  Initially  we  had 
planned  to  perform  the  jump  training  as  a  loading  model.  However,  if  jump  training  does  not 
show  a  significant  amount  of  bone  formation,  then  we  will  be  using  four-point  bending  as 
alternative  loading  regimen  and  real  time  PCR  as  the  end  point  measurement  of  the  bone 
formation  rate  in  two  strains  that  show  extreme  difference  in  bone  formation. 

This  report  includes  our  progress  for  the  last  twelve  months  of  our  proposal  (May  2002- 
May  2003). 
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Specific  Objective  1:  Optimize  a  method  that  can  induce  mechanical  strain  on  bone. 

Initially,  6-week  old  B6  mice  were  given  4-days  jump-up  and  jump-down  training  to 
evaluate  the  bone  formation  response.  The  mice  were  sacrificed  after  24hrs  of  the  last  training 
regimen  and  the  bone  formation  was  evaluated  in  terms  of  mRNA  expression  using  Real  time 
PCR.  In  this  experiment,  we  used  type-I  collagen  since  it  constitutes  approximately  95%  of  bone 
matrix  and  will  therefore  show  increased  expression  in  response  to  loading.  However,  we  found 
that  the  jump  training  did  not  produce  a  sufficient  amount  of  stress  on  the  bone  to  see  a 
significant  change  in  the  expression  of  type-I  collagen  between  control  and  experiment  (Table-1). 
The  values  in  the  table  represent  the  cycles  (Ct-values),  which  correspond  to  the  amount  of 
amplified  product  of  the  target  gene. 

Table-1  Quantification  of  type-1  collagen  expression  in  6- week  old  B6  mice  in  response 


Experiments 

Groups 

Cycles  ±  SD 

Gene 

0.1  Ug  of  RNA 

Fold  Change 

JUMP  DOWN-20cm 

Col  Al 

19.56  ±  1.03 

No  change 

Control 

Jump-down  20  cm 

Col  Al 

24.40  ±  0.97 

No  change 

Control 

P-actin 

23.49  ±  1.07 

No  change 

20cm  Jump-down 

P-actin 

22.69  ±  0.96 

No  change 

JUMP-DO  WN-40cm 

Control 

Col  Al 

19.83  ±0.72 

No  change 

Jump-down  40cm 

Col  Al 

19.40  ±0.7 

No  change 

Control 

P-actin 

21.98  ±  0.51 

No  change 

40cm  Jump-down 

P-actin 

22.75  ±  1.12 

No  change 

JUMP-UP  20cm 

Control 

Col  Al 

28.32  ±2.06 

Jump-up  20cm 

Col  Al 

28.65  ±  1.22 

No  change 

Control 

P-actin 

34.72  ±2.56 

No  change 

Jump-up  20cm 

P-actin 

35.51  ±  1.33 

No  change 

to  jump  training  using  real  time  PCR.  The  value  represents  Mean  ±  SD. 


N=4,  p>.  05  for  all  co-variance 
ColAl-type-1  collagen,  f-actin-beta-actin 

Therefore,  we  used  our  next  candidate  m-v/vo-loading  regimen,  namely  four-point 
bending.  The  load,  cycle,  and  frequency  used  for  the  study  were  similar  to  that  used  by  Reckers, 
et  al.,  (1998),  except  the  duration  of  training  was  changed  from  12  to  4-days  of  training  in  order 
to  evaluate  the  effect  of  mechanical  strain  on  bone  formation  in  a  short  duration.  The  loading 
regimen  consists  of  a  9N  load  at  2Hz  for  36cycles.  The  training  is  preformed  once  per  day.  The 
right  tibia  of  the  mouse  is  loaded  while  the  left  tibia  is  used  as  an  internal  control.  The  mice  are 
sacrificed  24  hrs  after  the  last  training  regimen  and  the  mRNA  message  is  quantified  using  Real 
time  PCR.  In  this  experiment,  in  order  to  evaluate  the  bone  formation  in  response  to  loading,  we 
used  two  major  bone  markers,  namely  type-I  collagen  and  osteocalcin.  In  addition  to  this,  we 
also  used  other  bone  formation  genes  namely  bone  sialoprotein,  alkaline  phosphate,  osteocalcin 
and  osteopontin. 

Four  days  four-point  bending  in  16-week  old  B6  and  C3H  mice  showed  a  significant 
(p<0.05)  increase  in  the  expression  of  bone  formation  genes  in  loaded  bones  as  compared  to 
unloaded  bones.  As  shown  in  Table-2&3,  mechanical  loading  resulted  in  two-fold  increase  in 
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type-I  collagen,  bone  sialoprotein,  and  alkaline  phosphates  in  both  B6  and  C3H  mice.  However, 
B6  mice  showed  1.4-fold  increase  in  Osteocalcin  while  C3H  did  not  show  any  change  between 
unloaded  and  loaded  bones. 

Table-2  Quantitative  Real  time  PCR  results  in  response  to  4-days  training  of  four-point 
bending  in  16- week  old  B6  mice  using  200  ng  RNA  concentration. 


Genes 

Groups 

Cycles  ±  SD 

P-value 

Fold  Change 

Type-I  Collagen 

Experiment 

20.17  ±1.37 

0.02 a 

2.2 

Control 

18.45  ±0.49 

Osteocalcin 

Control 

24.60  ±  0.86 

0.05 a 

1.42 

Experiment 

23.49  ±  0.70 

Bone  sialoprotein 

Control 

25.11  ±0.97 

0.01a 

2.2 

Experiment 

23.20  ±0.59 

Alkaline 

Phosphatase 

Control 

29.53  ±  1.34 

0.02a 

2.5 

Experiment 

27.60  ±  0.64 

Osteopontin 

Control 

24.01  ±  1.38 

0.38 

No  change 

Experiment 

23.36  ±0.77 

Actin 

Control 

19.19  ±1.04 

0.28 

Experiment 

18.60  ±0.49 

N  =  5„ap<0.05 

Table-3  Quantitative  Real  time  PCR  results  in  response  to  4-days  training  of  four-point 
bending  in  16- week  old  C3H  mice  using  200  ng  RNA  concentration. 


Genes 

Groups 

Cycles  ±  SD 

P-value 

Fold  Change 

Type-I  collagen 

Control 

23.77  ±  1.05 

0.0002* 

3.0 

Experiment 

20.47  ±0.58 

Osteocalcin 

Control 

26.16  ±0.86 

0.10 

No  change 

Experiment 

25.30  ±0.76 

Bone  sialoprotein 

Control 

27.06  ±  1.85 

0.02 a 

1.52 

Experiment 

24.74  ±  0.57 

Alkaline  Phosphatase 

Control 

31.10  ±  1.31 

0.003* 

2.0 

Experiment 

28.37  ±0.65 

Bone  Osteopontin 

Control 

26.92  ±  1.62 

0.25 

No  change 

Experiment 

25.89  ±0.93 

Actin 

Control 

22.82  ±  1.68 

0.07 

Experiment 

21.11  ±0.84 

N=  4,  ap<0.05,  *p<0.01 


Specific  Objective  2:  Determine  an  optimal  method  to  evaluate  the  bone  formation 

To  quantify  the  peak  gene  expression  in  response  to  loading,  we  used  a  technique  called 
real  time  PCR.  After  the  training  regimen,  mice  are  sacrificed  and  the  total  RNA  is  extracted  from 
the  loaded  and  unloaded  bone  by  Trizol  method.  This  RNA  is  used  to  synthesis  single  strand 
cDNA  by  reverse  transcription  assay  followed  by  Real  time  PCR  using  applied  biosystems.  The 
principle  of  this  Real  time  PCR  is  similar  to  an  ordinary  PCR  except  with  a  fluorescent  dye  which 
binds  to  amplified  product  and  emits  fluorescence  signals  that  are  viewed  in  terms  of  cycles  or 
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Ct-values.  The  Ct-values  correspond  to  the  amount  of  amplified  products  and  used  for  calculating 
the  mRNA  expression  of  the  target  gene.  In  order  to  validate  these  Ct-values,  we  also  run 
Universal  RNA  (Stratagene,  USA)  as  a  standard  during  each  reaction.  We  used  SYBR  Master 
mix  (Fluorescent  PCR  reagents)  kit  to  perform  the  above  reaction  in  applied  biosystems.  In  this 
four-point  bending  experiment,  to  evaluate  the  bone  formation  response  to  loading  we  used  two 
major  bone  markers,  namely  type-I  collagen  and  osteocalcin.  In  addition  to  this,  we  also  used 
other  bone  formation  genes,  namely  bone  sialoprotein,  alkaline  phosphatase,  and  osteopontin, 
since  they  play  a  significant  role  in  the  mineralization  of  bone.  B-actin  was  used  as  an  internal 
control  to  normalize  the  data  and  calculate  the  exact  fold  change  of  genes  by  applying  a 
mathematical  formula  (2'aaCT)  of  applied  biosystems. 

In  order  to  further  validate  the  four-point  bending  method  and  the  real  time  PCR  results, 
we  used  pQCT  and  histological  methods  that  are  well-established  endpoints  for  measuring  the 
bone  formation.  We  performed  12-days  training  of  four-point  bending  with  similar  magnitude  of 
load,  frequency  and  cycle  number  in  16-week  B6  and  C3H  mice.  Since  longer  periods  of  time  are 
required  to  see  new  bone  deposit,  we  used  12  days  of  training  instead  of  4  days.  As  shown  in 
Fig.l,  12-days  training  in  16-week  B6  mice  showed  maximum  change  in  total  area,  total  content, 
and  cortical  content  in  the  loaded  bones  compared  to  unloaded  bones.  Significant  changes  were 
also  observed  in  periosteal,  endosteal  circumference,  and  bone  density. 


i  UnLoaded  ■  Loaded 


1.6 


4  5  6  7  8  9 


Slice  number 


Fig.l  [a] 


i  UnLoaded  ■  Loaded 


2.5  - 


4  5  6  7  8  99 


Slice  number 


Fig.l  [b] 


Fig.l  [c] 


Fig.l  [d] 
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Fig.l  [e]  Fig.l  [f] 


Fig.l  Changes  in  bone  formation  parameters  in  response  to  12  days  training  of  four-point  bending 
measured  by  pQCT  in  tibia  of  16-week  female  B6  mice.  The  data  for  slice  4-9  which  represents  loading 
zone  are  shown  as  mean  ^standard  deviation  (n=9).  P-values  are  calculated  using  t-test  comparing  the 
loaded  and  unloaded  bones. 

(a)  Total  content  shows  pO.OOT* 

(b)  Total  area  shows  pO.OOl**  for  all  the  slices  except  slice-4  with  p<0.0T 

(c)  Periosteal  circumference  shows  p<0.001  for  all  the  slices  except  for  slice-4  with  pO.Ol* 

(d)  Endosteal  circumference  shows  pO.OOl**  for  slice  6,8,9  and  for  slice  5,  7  pO.Ol*  while  slice  4 
p>0.05 

(e)  Cortical  content  shows  p<0.001  for  all  the  slices 

(f)  Total  Density  showed  p<0.0T  for  all  the  slices 

Similar  to  the  B6,  C3H  mice  also  showed  maximum  change  in  total  area,  total  content, 
periosteal  and  endosteal  circumference  of  bone  in  response  loading.  However,  in  contrast  to  the 
B6,  there  was  no  change  in  total  density  in  the  loaded  bones  compared  to  unloaded  bone  as  shown 
in  Fig.2. 


Fig.2  [a]  Fig.2  [b] 
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*  UnLoaded 
■  Loaded 


4  5  6  7  8  9 

Slice  number 


Fig.2  [c] 


4  5  6  7  8  9 

Slice  number 


Fig.  2  [d] 
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Fig.2  [e]  Fig.2  [f] 

Fig.2  Changes  in  bone  formation  parameters  in  response  to  12  days  of  four-point  bending  measured  by 
pQCT  in  tibia  of  16-week  female  C3H  mice.  The  data  for  slice  4-9  which  represents  mean  +_standard 
deviation  (n=9).  P-values  are  calculated  using  t-test  comparing  the  loaded  and  unloaded  bones. 

a)  Total  Content  shows  p  <0.00 1  ** 

b)  Total  area  show  p<0.001**  for  all  the  slices  except  for  slice-4  p<0.05a 

c)  Periosteal  circumference  shows  pO.OOl**  except  for  slice-4  p<0.0T 

d)  Endosteal  circumference  shows  p<0.001**  for  slices  7, 9  while  slices  5,  8  showed  pO.Ol*.  Slices  4 
&6  showed  p>0.05 

e)  Total  bone  density  shows  p>0.05  for  all  the  slices 

f)  Cortical  Content  shows  pO.OOl**  for  all  the  slices  except  slice-4  with  p<0.01* 

Furthermore,  to  correlate  the  pQCT  and  gene  expression  results  using  real  time  PCR  in 
response  to  mechanical  loading,  we  performed  12  days  of  four-point  bending  in  retired  breeders 
of  B6  and  C3H  mice.  Twelve  days  training  of  four-point  bending  in  B6  retired  breeders  showed 
9-fold  increase  in  type-I  collagen  (ColAl)  and  7-fold  increase  in  alkaline  phosphatase  (ALP)  as 
compared  to  C3H  retired  breeders.  However,  the  bone  sialoprotein  showed  a  decrease  in  fold 
change  for  B6  compared  to  C3H.  This  suggests  that  long  durations  of  bending  treatment 
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decreased  the  expression  of  bone  sialoprotein  in  B6  (Fig.3).  Based  on  these  results,  type-I 
collagen  found  to  be  a  valuable  marker  that  is  increased  in  all  age  groups  in  response  to  bending. 
On  the  other  hand,  1 2  days  of  four-point  bending  using  pQCT  showed  increase  in  total  area,  total 
content,  periosteal,  and  endosteal  circumference  of  bone  for  B6  mice  compared  to  C3H  mice  of 
retired  breeders  as  shown  in  Fig.3a.  However,  for  the  total  density,  B6  showed  very  less  change 
while  C3H  showed  no  significant  change. 


*  B6  HC3H 


10  n  * 


ColAI  BSP  ALP 

Genes 


N=5,  *  p<0.01,  a p<0. 05 

Fig.3  Fold  change  in  expression  of  various  genes  in  response  to  12  days  of  four-point  bending  in  B6  and 
C3H  retired  breeders. 


"’pO-OOOl,  **p<0.001,  +p<0.05 


Fig.3a  Increase  in  bone  formation  parameters  in  response  to  12  days  training  of  four-point  bending 
measured  by  pQCT  in  tibia  of  B6  and  C3H  retired  breeders.  The  data  shows  percentage  of  increase  in  TA- 
Total  Content,  TA-Total  Area,  PC-Periosteal  circumference,  EC-Endosteal  circumference  as  Mean  + 
Standard  deviation  of  loaded  zone.  P-values  are  calculated  by  comparing  the  percentage  of  new  bone 
formation  between  B6  and  C3H  mice. 


We  used  Real-time  PCR  and  pQCT  as  two  methods  for  measuring  the  bone  formation 
rate  in  response  to  four-point  bending.  Fig.3b  shows  the  correlation  between  Bone  Mineral 
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Density  and  Type-I  Collagen  expression.  As  shown,  there  is  a  significant  correlation  (p<0.05) 
between  the  BMD  [Measured  by  pQCT]  and  collagen  expression  [Measured  by  real  time  PCR]. 
Based  on  these  results,  we  conclude  that  Real  time  PCR  is  the  most  sensitive  and  fastest  method 
to  evaluate  the  bone  formation  in  response  to  four-point  bending.  Thus,  we  will  be  using  Real 
Time  PCR  as  the  endpoint  in  our  QTL  study. 


Fig.3b  Correlation  between  bone  mineral  density  and  type-I  collagen  expression  in  response  to  four-point 
bending  in  different  age  groups  of  B6  mice.  The  x-axis  represents  BMD  change  measured  by  percentage 
increase  in  loaded  bone  compared  to  corresponding  contralateral  bone.  The  y-axis  shows  change  in 
collagen  gene  expression  in  response  to  four-point  bending  as  measured  by  MCt  for  loaded  bone  vs 
unloaded  bone  using  b-actin  as  control  gene. 

Specific  Objective  3:  Determine  an  optimal  duration  of  Training  period 

Ten-week  old  B6  mice  were  subjected  to  2  days,  4  days  and  8  days  training  in  order  to 
select  an  optimal  duration  of  training  that  show  maximum  effect  on  bone  formation.  2  days  of 
four-point  bending  in  B6  mice  did  not  produce  any  significant  effect  on  bone.  The  expression  of 
bone  formation  genes  in  loaded  bones  was  similar  to  the  control  shown  in  Table-4. 

Table  -4  Quantitative  Real  time  PCR  results  in  response  to  2  days  training  of  four-point 
bending  in  1 0-week  old  B6  mice  using  200  ng  RNA  concentration. 


Genes 

Groups 

Cycles  ±  SD 

P-value 

Fold  Change 

Type-I  collagen 

Control 

20.61  ±1.9 

0.97 

No  change 

Experiment 

20.58  ±  0.4 

Osteocalcin 

Control 

24.85  ±  2.2 

0.52 

No  change 

Experiment 

24.17  ±0.6 

Bone  sialoprotein 

Control 

26.59  ±  1.4 

0.82 

No  change 

Experiment 

26.76  ±  0.59 

Alkaline  Phosphatase 

Control 

30.70  ±  1.9 

0.86 

No  change 

Experiment 

30.86  ±0.5 

Actin 

Control 

22.40  ±  2.4 

0.48 

Internal  control 

Experiment 

21.60  ±0.49 

N=5,  p>0.05 
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Four  days  four-point  bending  in  10-week  old  B6  mice  resulted  in  2.9  and  4.5  fold 
increase  in  type-1-  collagen  and  bone  sialoprotein  respectively.  On  the  other  hand,  no  change  was 
observed  in  alkaline  phosphates,  osteocalcin  and  osteopontin  (Table-5).  B-actin  was  used  as 
internal  control  to  normalize  the  data  in  order  to  obtain  the  fold  change  of  the  target  gene. 

Table-5  Quantitative  Real  time  PCR  results  in  response  to  4  days  four-point  bending  in 
10-week  old  B6  mice  using  200  ng  RNA  concentration. 


Genes 

Groups 

Cycles  ±  SD 

P-value 

Fold  Change 

Type-I  Collagen 

Experiment 

21.75  ±  1.8 

0.06 

2.9 

Control 

19.45  ±  1.5 

Bone  sialoprotein 

Control 

27.21  ±  1.4 

0.02a 

4.5 

Experiment 

24.22  ±1.9 

Osteocalcin 

Control 

25.28  ±  1.1 

0.49 

No  change 

Experiment 

24.76  ±1.06 

Alkaline  phosphatase 

Control 

31.88  ±2.4 

0.14 

No  change 

Experiment 

29.55  ±2.1 

Osteopontin 

Control 

24.36  ±  1.3 

0.34 

No  change 

Experiment 

23.57  ±  1.1 

Actin 

Control 

22.56  ±  1.6 

0.47 

Internal  Control 

Experiment 

21.84  ±  1.3 

N=5,  a  p<  0.0 5 


Eight  days  of  four-point  bending  in  tibia  of  10-week  old  B6  mice  showed  significant 
change  in  expression  of  bone  formation  genes  compared  to  2-days  and  4-days  of  training.  As 
shown  in  Table-6,  8-days  of  training  showed  4-fold  increase  in  type-I  collagen,  2.7-fold  in 
osteocalcin,  3.6  fold  in  Bone  sialoprotein  and  2.8  fold  in  alkaline  phosphotase. 


Table-6  Quantitative  Real  time  PCR  results  in  response  to  8  days  four-point  bending  in 
10- week  old  B6  mice  using  200  ng  RNA  concentration. 


Genes 

Groups 

Cycles  ±  SD 

P-value 

Fold  Change 

Type-I  collagen 

Control 

22.02  ±  0.4 

3.8 

Experiment 

19.35  ±0.4 

Osteocalcin 

Control 

24.94  ±1.1 

0.01a 

2.7 

Experiment 

22.76  ±  0.9 

Bone  sialoprotein 

Control 

26.79  ±0.8 

0.003’ 

3.7 

Experiment 

24.16  ±  1.1 

Alkaline  Phosphatase 

Control 

29.78  ±  0.9 

0.001* 

2.9 

27.53  ±  0.5 

Actin 

22.56  ±  1. 1 

0.24 

Internal  control 

Experiment 

21.83  ±0.6 

N=5,  *p<0.01,  ** p<0.001 ,  ap<0.05 


12 


Specific  Objective  4:  Determine  optimal  mouse  strain  pair  that  show  extreme 
difference  in  bone  formation. 

With  the  optimized  loading  regimen,  four-point  bending  was  carried  out  in  1 0-week-old 
B6,  Balbc,  AKR/J,  129J,  NZB  mice  and  the  bone  formation  was  evaluated  using  Real  time  PCR. 
AKR/J  showed  one-fold  change  in  type-I  collagen,  alkaline  phosphatase,  and  two-fold  change  in 
Bone  sialoprotein  while  no  change  in  osteocalcin  was  seen  (Table-7).  129J  showed  3  fold  change 
in  type-I  collagen  and  bone  sialoprotein.  However,  no  change  in  alkaline  phosphatase, 
osteocalcin,  and  osteopontin  was  observed  (Table-8).  NZB  showed  no  change  in  any  of  the  bone 
markers  in  response  to  four-point  bending  (Table-9).  In  the  case  of  Balbc,  the  bone  (tibia)  broke 
at  9N  load,  therefore,  a  6N  load  was  applied.  There  was  no  change  in  expression  of  bone 
formation  genes  between  control  and  experiment  in  Balbc  mice  (Table- 10).  The  B6  mice  showed 
three-four  fold  increase  in  type-I  collagen  and  bone  sialoprotein.  However,  there  was  no  change 
in  alkaline  phosphatase  and  osteocalcin  (Table-11).  Overall,  AKR,  129J  and  B6  mice  showed 
similar  gene  expression  in  response  to  four-point  bending  while  Balbc  and  NZB  showed  no 
significant  change  in  expression  of  bone  formation  genes  shown  in  Table-12 


Table-7  Quantitative  Real  time  PCR  results  in  response  to  4-days  four-point  bending  in 
10-week  old  AKR/J  using  200  ng  RNA  concentration. 


Genes 

Groups 

Cycles  ±  SD 

P-value 

Fold  Change 

Type-I  Collagen 

Experiment 

18.16  ±0.9 

0.01 a 

1.0 

Control 

19.35  ±0.4 

Bone  sialoprotein 

Control 

25.27  ±  1.7 

0.01a 

2.7 

Experiment 

22.76  ±1.1 

Alkaline  Phosphatase 

Control 

24.23  ±0.8 

0.02 a 

1.3 

Experiment 

22.87  ±  0.9 

Osteocalcin 

Control 

22.11  ±0.4 

0.44 

No  change 

Experiment 

21.74  ±  1.0 

Actin 

Control 

20.73  ±  1.5 

0.1 

Internal  control 

Experiment 

19.65  ±0.9 

N=5,ap<0.05 

Table-8  Quantitative  Real  time  PCR  results  in  response  to  4-days  four-point  bending  in 
10-week  old  129J  using  200  ng  RNA  concentration. 
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Genes 

Groups 

Cycles  ±  SD 

P-value 

Fold  Change 

Type-I  Collagen 

Experiment 

20.87  ±0.7 

0.06 

2.6 

Control 

19.66  ±0.9 

Bone  sialoprotein 

Control 

23.48  ±  0.3 

0.05 a 

2.8 

Experiment 

22.14  ±  1.3 

Alkaline  Phosphatase 

Control 

27.33  ±  1.3 

0.57 

No  change 

Experiment 

26.80  ±  1.5 

Osteopontin 

Control 

29.33  ±  0.4 

0.37 

No  change 

Experiment 

29.93  ±  1.3 

Osteocalcin 

Control 

26.49  ±  0.6 

0.53 

No  change 

Experiment 

26.26  ±  0.4 

Actin 

Control 

19.62  ±1.1 

0.53 

Internal  Control 

Experiment 

20.07  ±  1.0 

N=5,  ap<0.05 


Table-9  Quantitative  Real  time  PCR  results  in  response  to  4-days  four-point  bending  in 
10-week  old  NZB  using  200  ng  RNA  concentration. 


Genes 

Groups 

Cycles  ±  SD 

P-value 

Fold  Change 

Type-I  Collagen 

Experiment 

19.50  ±2.4 

0.63 

No  change 

Control 

18.94  ±  1.3 

Bone  sialoprotein 

Control 

23.89  ±  1.6 

0.12 

No  change 

Experiment 

24.77  ±2.6 

Osteocalcin 

Control 

22.16  ±  1.7 

0.50 

No  change 

Experiment 

24.13  ±2.3 

Actin 

Control 

20.86  ±  1.4 

0.31 

Internal  control 

Experiment 

21.92  ±  1.9 

N=5,  p>0.05 


Table-10  Quantitative  Real  time  PCR  results  in  response  to  4-days  four-point  bending  in 
BalbclO-week  old  mice  using  200  ng  RNA  concentration. 


Genes 

Groups 

Cycles  ±  SD 

P-value 

Fold  change 

Type-I  Collagen 

Experiment 

17.80  ±  1.4 

0.98 

No  change 

Control 

1 7.78  ±  1.7 

Bone  sialoprotein 

Control 

23.32  ±  1.5 

0.49 

No  change 

Experiment 

22.59  ±  1.5 

Alkaline  Phosphatase 

Control 

25.32  ±2.1 

0.81 

No  change 

Experiment 

25.00  ±  1.4 

Osteopontin 

Control 

28.51  ±2.0 

0.88 

No  change 

Experiment 

28.35  ±  1.1 

Actin 

Control 

17.86  ±  1.0 

0.33 

Internal  control 

17.21  ±  0.9 
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N=5p>0.05 


Table-11  Quantitative  Real  time  PCR  results  in  response  to  4  days  training  of  four-point 
bending  in  10- week  old  B6  mice  using  200ng  RNA  concentration 


Genes 

Groups 

Cycles  ±  SD 

P-value 

Fold  Change 

Type-I  Collagen 

Experiment 

21.75  ±  1.8 

0.06 

2.9 

Control 

19.45  ±  1.5 

Bone  sialoprotein 

Control 

27.21  ±  1.4 

0.02  a 

4.5 

Experiment 

24.22  ±  1.9 

Osteocalcin 

Control 

25.28  ±  1.1 

0.49 

No  change 

Experiment 

24.76  ±  1.0 

Alkaline  phosphatase 

Control 

31.88  ±  2.4 

0.14 

No  change 

Experiment 

29.55  ±2.1 

Osteopontin 

Control 

24.36  ±  1.3 

0.34 

No  change 

Experiment 

23.57  ±  1.1 

Actin 

Control 

22.56  ±  1.6 

0.47 

Internal  control 

Experiment 

21.84  ±  1.3 

N=5,  ap-<0. 05 


Table-12  Ranking  of  mouse  strain  [Good  and  Poor  responder]  based  upon  the  expression 
of  bone  formation  genes  in  response  to  4  days  of  four-point  bending. 


Bone  Formation 
Markers 

Fold  change  in  response  to  4-days  training  of  four-point 
bending  in  different  inbred  strains  of  mice  [10-week  old] 

Good  Responder  1 

|  Poor  Responder 

Genes 

_B6 

129J 

AKR/J 

NZB 

Balbc 

Type-I  collagen 

2.9 

2.6 

1.0 

. 

Bone  sialoprotein 

4.5 

2.8 

2.7 

. 

Alkaline  phosphatase 

1.3 

Osteocalsin 

_ 

_ 

_ 

_ 

- 

-No  change 


Specific  Objective  5:  Determine  appropriate  age  of  mice  that  shows  significant 
affect  on  bone  formation  through  optimized  protocol. 

With  the  optimized  loading  regimen  [9N  load  at  2Hz  -36cycles]  four-point  bending  was 
applied  B6  mice  of  10-week  old,  16-week  old  and  Retired  Breeders  and  the  bone  formation  was 
evaluated  using  real  time  PCR.  As  shown  in  Table-13,  10-week  old  mice  in  response  to  4-days 
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training  showed  increased  expression  of  type-I  collagen  and  Bone  sialoprotein  (p-<0.05)  while  no 
change  in  osteocalcin  and  alkaline  phosphates.  While  1 6-week  old  mice  showed  two-fold  change 
in  type-I  collagen,  bone  sialoprotein,  alkaline  phosphates  and  1.4  fold  change  in  osteocalcin 
shown  in  Table-14.  Interestingly  in  the  retired  breeders  we  found  that  increase  in  age  show  a 
better  response  of  bone  formation  in  response  to  4-days  training  of  9N  load.  The  type-I  collagen 
and  Bone  sialoprotein  showed  a  significant  fold  change  with  p<0.01  while  alkaline  phosphotase 
and  osteocalcin  showed  p<0.05  (Table- 15). 

Table-13  Quantitative  Real  time  PCR  results  in  response  to  4-days  four-point  bending  in 
10-week  old  female  B6  mice  using  200  ng  RNA  concentration 


Genes 

Groups 

Cycles  ±  SD 

P-value 

Fold  Change 

Type-I  Collagen 

Experiment 

21.75  ±  1.8 

0.06 

2.9 

Control 

19.45  ±  1.5 

Bone  sialoprotein 

Control 

27.21  ±  1.4 

0.02 a 

4.5 

Experiment 

24.22  ±  1.9 

Osteocalcin 

Control 

25.28  ±  1.1 

0.49 

No  change 

Experiment 

24.76  ±  1.06 

Alkaline  phosphatase 

Control 

31.88  ±2.4 

0.14 

No  change 

Experiment 

29.55  ±2.1 

Osteopontin 

Control 

24.36  ±  1.3 

0.34 

No  change 

Experiment 

23.57  ±  1.1 

Actin 

Control 

22.56  ±  1.6 

0.47 

Internal  control 

Experiment 

21.84  ±  1.3 

N=5,  ap<0. 05 


Table-14  Quantitative  Real  time  PCR  results  in  response  to  4-days  four-point  bending  in 
1 6-week  old  female  B6  mice  using  200  ng  RNA  concentration. 


Genes 

Groups 

Cycles  ±  SD 

P-value 

Fold  Change 

Type-I  Collagen 

Experiment 

20.17  ±  1.37 

0.02 a 

2.2 

Control 

18.45  ±0.49 

Osteocalcin 

Control 

24.60  ±0.86 

0.05 a 

1.42 

Experiment 

23.49  ±0.70 

Bone  sialoprotein 

Control 

25.11  ±0.97 

0.01 a 

2.2 

Experiment 

23.20  ±0.59 

Alkaline 

Phosphatase 

Control 

29.53  ±  1.34 

0.02 a 

2.5 

Experiment 

27.60  ±  0.64 

Osteopontin 

Control 

24.01  ±1.38 

0.38 

No  change 

Experiment 

23.36  ±0.77 

Actin 

Control 

19.19  ±  1.04 

0.28 

Internal  control 

Experiment 

18.60  ±0.49 

N=  5,  dp<0.5 
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Table-15  Quantitative  Real  time  PCR  results  in  response  to  4  days  four-point  bending  in 
B6  Retired  breeders  using  200  ng  RNA  concentration 


Genes 

Groups 

Cycles  ±  SD 

P-value 

Fold  Change 

Type-I  collagen 

22.26  ±  1.0 

0.0007’ 

6.5 

18.61  ±0.54 

Bone  sialoprotein 

Ej 

25.58  ±  1.1 

0.002* 

6.5 

21.93  ±0.98 

Alkaline  phosphotase 

Control 

30.16  ±  1.5 

0.02 a 

4 

Experiment 

27.14  ±1.2 

Osteocalcin 

Control 

27.28  ±  1.0 

0.04 a 

2 

Experiment 

25.33  ±  1.1 

Actin 

Control 

20.82  ±  0.64 

0.12 

Internal  control 

Experiment 

19.88  ±0.82 

N  =  4,  ap-<0. 05,  *p-<0.01 


Key  Findings 

1.  Determination  of  an  appropriate  response  parameter.  Because  pQCT  is  a  time- 
consuming  labor-intensive  method,  we  sought  a  more  rapid,  but  yet  sufficiently  sensitive, 
method  for  measuring  bone  formation  and  response  to  mechanical  loading.  Measurement 
of  bone  marker  expression  by  real  time  PCR  has  not  been  previously  used  to  evaluate 
bone  formation  in  a  quantitative  manner.  Therefore,  it  was  important  for  us  to  correlate 
the  changes  within  an  acceptable  bone  formation  measurement,  namely  pQCT  with  our 
data  obtained  from  real  time  PCR  on  gene  expression.  In  this  study,  we  found  a  good 
correlation  between  the  real  time  PCR  changes  and  the  pQCT  changes  (pQCT  is  a 
measure  of  bone  density  whereas  we  have  interpreted  this  to  indicate  bone  formation 
simply  because  the  amount  of  change  in  bone  density  is  a  reflection  of  bone  formation). 
In  conclusion,  we  found  that  real  time  PCR  evaluating  specific  bone  marker  genes 
fulfilled  our  criteria  for  a  sensitive  reproducible  and  producible  measurement  that  shows 
a  nice  response  to  four-point  bending. 

2.  Duration  of  application  of  mechanical  strain.  Our  studies  indicate  that  four  days  of  four- 
point  bending  in  our  femoral  in  vivo  loading  regimen  is  adequate  for  a  measurable 
response  in  our  proposed  QTL  studies. 

3.  Identification  of  bone  response.  These  are  genes  that  can  be  assayed  as  a  group  in  real 
time  PCR  for  assessment  of  the  bone  formation  response  to  four-point  bending.  We 
found  that  alkaline  phosphatase,  osteocalsin,  bone  sialoprotein,  and  type-1  collagen  genes 
showed  statistically  significant  increases  in  RNA  harvested  from  bone  tissue  following 
four  days  of  four-point  bending.  Consequently,  these  genes  were  selected  for  further 
application. 

4.  Determination  of  the  animal  chronological  age  for  application  of  four-point  bending.  We 
found  that  retired  breeder  mice  and  1 6-week  old  mice  showed  the  greatest  increase  in  our 
bone  formation  parameters  in  response  to  four-point  bending.  However,  we  would  prefer 
to  use  mice  10  weeks  of  age  since  we  could  progress  in  our  QTL  studies  much  more 
rapidly  by  using  the  younger  mice.  We  are  now  in  the  process  of  settling  whether  we  can 
use  the  10-week  old  mice  for  our  four-point  bending  QTL  study. 
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5.  The  two  strains  of  mice  most  appropriate  for  the  four-point  bending  OTL  study.  In 
addition  to  C3  and  B6  animals  we  have  studied  four  additional  strains;  AKR,  129J,  NZB, 
and  the  Balb/c  mice.  The  AKR  and  the  129J  showed  a  good  response  whereas  the  NZB 
and  the  Balb/c  showed  no  response.  We  are  now  in  the  process  in  determining  whether 
some  combination  of  these  latter  four  strains  would  be  superior  to  the  C3  and  B6  for  the 
high  and  low  response  pair. 

Reportable  Outcomes 

None 

Conclusions 

1.  Goals  to  accomplish  before  beginning  the  OTL  studies.  We  have  now  completed  all  of 
the  goals  required  for  QTL  studies  except  one,  which  we  did  not,  include  in  our  technical 
objective  namely  a  dose  response.  Apart  from  that,  we  are  ready  to  begin  the  QTL  study. 
We  are  now  in  the  process  of  doing  a  dose  response  to  determine  the  optimal  response 
range  which  best  distinguishes  the  two  different  strains  of  mice,  the  C3H  mice  and  B6 
mice. 

2.  Selection  of  the  optimal  strains  of  mice.  Based  on  the  data  so  far,  the  C3H  and  B6  mice 
are  the  best  strain  pair  and  as  mentioned  above  the  final  decision  on  this  issue  will  be 
based  on  the  dose  response. 

3.  Real  time  PCR  measurements  of  bone  markers  as  a  surrogate  for  pOCT  measurement. 
When  we  plotted  type-I  collagen  as  our  real  time  PCR  measurement  and  the  increment  in 
pQCT  bone  density  over  the  unloaded  bone  control  bone,  we  found  a  significant 
correlation,  which  supports  the  view  that  we  can  use  real  time  PCR  as  a  much  faster 
means  to  obtain  our  phenotype  after  four-point  bending. 

4.  Additional  phenotypes.  Not  only  can  we  use  these  measurements  of  message  level  of 
bone  proteins  as  indices  of  the  bone  formation  response,  we  can  also  use  them  separately 
to  evaluate  their  level  of  expression  as  a  separate  phenotypes  to  determine  the  QTLs 
responsible  for  quantitative  level  of  expression  of  each  real  time  PCR  mRNAs. 
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2.  Molecular  Genetic  Studies  on  Bone  Mechanical  Strain  -  in  vitro  studies 
Introduction 

It  is  well  established  that  mechanical  loading  leads  to  an  increase  in  bone  density  and  that 
immobilization  leads  to  a  loss  in  bone  density.  Although  the  response  of  bone  cells  to  mechanical 
stimuli  is  relatively  well  understood,  the  knowledge  about  the  signaling  mechanisms  and  the 
genes  involved  in  the  mechanical  regulation  of  bone  structure  and  function  are  less  understood  or 
limited.  In  our  recent  study  we  have  shown  that  two  inbred  strains  of  mice  (C3H  and  B6)  showed 
differing  responses  in  bone  turnover  to  mechanical  loading,  suggesting  that  the  genetic 
component  of  bone  response  to  mechanical  loading  must  be  biologically  significant.  To  identify 
the  genetic  component  involved,  our  goal  is  to  apply  a  powerful  genetic  approach  known  as 
quantitative  trait  loci  (QTL)  in  combination  with  in  vitro  studies  on  the  biochemical  pathways 
involved  in  mechanical  stress  signaling  in  order  to  determine  the  genes  responsible  for  the 
mechanical  stress  differential  responses  between  the  two  inbred  strains  of  mice  that  show  good 
and  poor  response  of  bone  formation. 

Technical  Objectives: 

To  evaluate  the  effects  of  mechanical  signaling  using  a  physiologically  relevant 
CytoDyne  flow  chamber  to  produce  a  fluid  flow  shear  strain  for  evaluation  of  proliferation  and 
differentiation  and  also  for  studies  of  gene  expression  and  signal  transduction  pathways  in 
cultures  of  C3H  and  B6  mouse  osteoblasts.  We  will  have  the  following  specific  objectives  during 
the  second  year  of  this  grant  period: 

1)  We  will  complete  our  initial  evaluation  of  strain-induced  tyrosine  phosphorylation 
levels  of  key  signaling  proteins  by  immunoprecipitation-immunoblotting  approach 
using  available  antibodies.  This  technology  will  be  applied  to  compare  the  responses 
of  C3H  vs.  B6  mice. 

2)  Apply  in-house  microarrays  to  evaluate  changes  in  gene  expression  of  those  genes 
known  to  be  involved  in  mediating  loading  signal. 

3)  The  genes  that  we  identified  by  the  above  technique  (i.e.,  the  microarray  technique) 
will  be  further  compared  with  real  time  PCR. 

In  our  previous  report,  we  optimized  the  conditions  to  develop  an  in  vitro  system  to  apply  defined 
mechanical  strain  to  cultured  bone  cells  and  compared  the  phenotypic  differences  between  the 
C3H  and  B6  mice  in  response  to  fluid  flow  shear  stress.  We  found  that  osteoblasts  isolated  from 
B6  mice  were  highly  responsive  to  shear  stress  in  terms  of  cell  proliferation  and  differentiation. 
On  the  other  hand,  osteoblasts  isolated  from  C3H  mice  were  unresponsive  to  the  same  shear 
stress.  The  following  is  our  progress  toward  each  of  our  objectives  during  Year  2  of  this  period 
(i.e.,  the  last  12  months). 

Specific  Objective  1:  Evaluate  strain-induced  tyrosine  phosphorylation  levels  of  key 
signaling  proteins  in  C3H  and  B6  mice. 

In  order  to  clarify  the  signaling  pathways  which  could  lead  to  the  changes  in  the  fluid 
flow  induced  cell  proliferation  and  differentiation  in  B6  bone  cells  and  not  in  the  bone  cells 
isolated  from  C3H  mice,  we  studied  the  fluid  flow  induced  phosphorylation  of  MAPK  and 
integrin  expression  in  these  cells,  since  the  MAPK  pathway  and  the  integrin  pathway  have  each 
been  implicated  in  the  mechanical  signaling  mechanism.  As  shown  in  Figure  1,  fluid  flow  did  not 
have  any  significant  effect  on  the  integrin  pi  expression  in  the  C3H  bone  cells.  Interestingly,  B6 
bone  cells  showed  a  significant  increase  in  the  integrin  pi  expression  when  subjected  to  fluid 
flow.  These  findings  suggest  that  shear  stress-  induced  osteoblast  proliferation  and  differentiation 


in  B6  mice  may  be  associated  with  the  shear  stress-dependent  upregulation  of  integrin  pi 
expression. 
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Figure  1.  Effect  of  shear  stress  on  integrin  pi  in  osteoblasts  isolated  from  C3H  and  B6  mice.  Cells 
were  subjected  to  fluid  flow  of  20  dynes/cm2  for  30  minutes.  (A).  Cell  lysates  were  immunoblotted  with 
anti-integrin  pi  and  anti-actin  antibodies.  (B).  The  graph  represents  the  densitometric  measurements  of 
integrin  pi  levels  from  western  blots  normalized  by  actin. 

Further,  when  the  effect  of  fluid  flow  on  MAPK  phosphorylation  was  studied,  it  was 
observed  that  there  was  a  significant  increase  in  the  phosphorylation  levels  of  both  the  ERK1  and 
ERK2  in  B6  bone  cells  when  subjected  to  fluid  flow.  However,  in  C3H  bone  cells,  no  change 
was  noted  either  in  ERK1  or  ERK2  phosphorylation  levels  in  response  to  fluid  flow  (Figure  2). 
These  findings  suggest  that  the  ERK  signaling  pathway  is  essential  in  the  osteogenic  response  to 
mechanical  stimuli  in  B6  bone  cells.  However,  a  remaining  key  question  is  whether  activation  of 
either  ERK1  or  ERK2  or  both  is  required. 


Figure  2.  Effect  of  shear  stress  on  phosphorylation  levels  of  MAPK  in  osteoblasts 
isolated  from  C3H  and  B6  mice.  Cells  were  subjected  to  fluid  flow  of  20  dynes/cm2  for  30  minutes.  (A): 
Cell  lysates  were  immunoblotted  with  anti-phospho-ERK  and  anti-actin  antibodies.  (B&C):  The  graph 
represents  the  densitometric  measurements  of  ERK1  and  ERK2  phosphorylation  levels  from  western  blots 
normalized  by  actin. 

To  further  evaluate  the  involvement  of  the  ERK-signaling  pathway  as  to  whether  the 
protein  kinase  activity  of  ERK1  and/or  2  is  essential  for  the  fluid  flow  shear  stress-induced 
osteoblast  proliferation,  we  overexpressd  dominant-negative  constructs  of  ERK1  and  ERK2  in 


bone  cells,  and  determined  whether  blocking  the  activation  of  the  ERK1  and/or  ERK2  would 
prevent  the  mechanical  stress-induced  bone  cell  proliferation. 

Shear  stress-induced  cell  proliferation  was  studied  in  human  TE85  osteosarcoma  cells 
that  were  transduced  with  MLV  retroviral-based  vectors  expressing  either  the  wild  type  (wt)  or 
dominant  negative,  kinase-dead  (kd)  ERK1  and  2.  Due  to  very  low  transduction  efficiency  in 
osteoblasts  isolated  from  B6  mice,  we  used  human  TE85  cells  in  this  part  of  the  study  as  the 
transduction  efficiency  of  TE85  cells  is  very  high  (79-90%).  A  HA  tag  was  added  to  the  N- 
terminus  of  each  construct  to  help  to  distinguish  the  overexpressed  enzyme  from  the  endogenous 
enzyme.  In  order  to  distinguish  the  overexpressed  ERK2  from  the  endogenous  protein,  a  tag  of  35 
amino  acids  in  size  was  introduced  to  the  C-terminus  of  the  protein.  An  MLV-red  fluorescent 
protein  (REP)  vector  was  also  included  as  a  control  for  comparison. 

To  ascertain  that  the  wild-type  ERK  expression  vectors  expressed  functionally  active 
ERK1  and  ERK2,  and  that  the  kinase-dead  mutants  of  ERK  vectors  expressed  inactive  ERK1  and 
ERK2,  protein  kinase  activity  was  assayed.  As  shown  in  Figure  3,  ERK1  and  ERK2  expressed  by 
each  corresponding  wild-type  expression  vector  (El  and  E2)  were  active  as  a  MAPK  whereas 
ERK1  and  ERK2  expressed  by  each  corresponding  kinase-dead  expression  vector  (el  and  e2) 
were  completely  inactive  as  a  MAPK.  Furthermore,  the  results  here  also  suggest  that  the  addition 
of  a  35  amino  acid  tag  at  the  C-terminus  did  not  appear  to  have  any  effect  on  the  protein  kinase 
activity  of  ERK2. 


TE85  El  E2  el  e2  Neg 

Control 


Figure  3.  The  Protein  kinase  activity  of  protein  products  expressed  by  each  ERK  expression  vector. 
(A):To  ascertain  that  the  wild-type  ERK  expression  vectors  expressed  functionally  active  ERK1  and 
ERK2,  respectively,  and  that  the  kinase-dead  mutants  of  ERK  vectors  expressed  inactive  ERK1  and  ERK2, 
respectively,  the  corresponding  ERK  in  lysates  of  E.  coli  expressing  wild  type  GST-ERK1  and  GST-ERK2 
vectors  (El  and  E2,  respectively)  or  that  of  E.  coli  expressing  kinase-dead  GST-ERK1  and  ERK2  vectors 
(el  and  e2,  respectively)  were  isolated  with  glutathione-conjugated  beads.  The  protein  kinase  assay  was 
then  assayed  with  MBP  as  the  substrate.  A  negative  control  (lysates  of  E.  coli  without  the  ERK  construct) 
and  TE85  cell  lysate  (positive  control)  were  included  for  comparison.  (B):The  graph  represents  the 
densitometric  measurements  of  kinase  activity. 

Human  TE85  osteosarcoma  cells  were  transduced  three  times  with  each  of  the  test 
retroviral  vectors  with  a  total  MOI  of  30.  To  ensure  that  human  TE85  cells  transduced  with  ERK 
expressing  MLV  vectors  indeed  overexpress  corresponding  ERKs,  we  measured  the  protein  level 
of  overexpressed  (oERKl  and  oERK2,  respectively)  and  endogeneous  (eERKl  and  eERK2) 
ERK1  and  2  in  cell  lysates  of  transduced  TE85  cells.  As  shown  in  Figure  4,  TE85  cells 
transduced  with  MLV  vectors  expressing  either  wild  type  or  kinase-dead  mutants  of  ERK  1  or  2 
led  to  overexpression  of  ERK  1  and  ERK2  protein,  respectively. 


Figure  4.  Overexpression  of  ERK  proteins  in  human  TE85  cells  transduced  with  MLV  vectors 
expressing  wild-type  (El  and  E2)  or  kinase-dead  mutants  (el  and  e2)  of  ERK.  (A):To  ensure  that 
human  TE85  cells  transduced  with  ERK  expressing  MLV  vectors  indeed  overexpress  corresponding  ERKs, 
the  protein  level  of  overexpressed  (oERKl  and  oERK2,  respectively)  and  endogeneous  (eERKl  and 
eERK2)  ERK1  and  2  was  measured  in  cell  lysates  of  TE85  cells  transduced  with  respective  MLV  vectors 
by  western  blot  analysis.  The  oERK2  has  an  extra  35  aa  at  the  C-terminus  compared  to  eERK2.  Cells 
transfected  with  MLV-RFP  vector  were  included  as  a  negative  control.  The  blot  was  stripped  and  reblotted 
with  anti-actin  antibody  to  assess  the  protein  loading.(B):  The  graph  represents  the  densitometric 
measurements  of  ERK1  and  ERK2  expression  levels  from  western  blots  normalized  by  actin. 

To  test  whether  the  protein  kinase  activity  of  ERK1  and/or  ERK2  is  essential  for  the  fluid 
flow  shear  stress-induced  proliferation  of  osteoblasts,  we  determined  the  effect  of  overexpression 
of  kinase-dead  ERK1  (el)  or  ERK2  (e2)  on  [3H]thymidine  incorporation  in  response  to  shear 
stress.  Figure  5  shows  that  TE85  cells  exhibited  a  significant  (~2-fold)  increase  in  cell 
proliferation  in  response  to  a  30-min  steady  20  dynes/cm2  shear  stress.  Furthermore, 
overexpression  of  either  wild  type  ERK1  or  ERK2  resulted  in  a  similar  increase  in  cell 
proliferation  in  response  to  shear  stress  as  compared  to  the  controls.  On  the  other  hand,  either  a 
dominant-negative,  kinase-dead  ERK1  (el)  or  ERK2  (e2)  mutant  completely  blocked  the 
mitogenic  response  of  TE85  cells  to  the  shear  stress,  suggesting  that  ERK1  and  ERK2  each  have 
an  important  role  in  mediating  the  shear  stress-induced  bone  cell  proliferation.  Figure  5  also 
shows  that  MLV-RFP  showed  a  similar  response  to  shear  stress  as  normal  TE85  cells  thus 
indicating  that  transfection  of  TE85  cells  with  the  MLV  vector  did  not  affect  the  bone  cell 
mitogenic  response  of  these  cells  to  the  30-min  shear  stress. 
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Figure  5:  Effect  of  overexpression  of  ERK1  (El  and  el)  or  ERK2  (E2  and  e2)  on  the  shear  stress- 
induced  increase  in  |3H]thymidine  incorporation.  To  test  whether  the  protein  kinase  activity  of  ERK1 
and/or  ERK2  is  essential  for  the  fluid  flow  shear  stress-induced  proliferation  of  osteoblasts,  we  determined 
the  effect  of  overexpression  of  kinase-dead  ERK1  (el)  or  ERK2  (e2)  on  [3H]thymidine  incorporation  in 
response  to  shear  stress.  *  *p<0.01. 

These  results  clearly  show  that  the  protein  kinase  activity  of  both  ERK1  and  ERK2  is 
essential  for  the  fluid  flow  shear  stress-induced  osteoblast  proliferation.  The  most  intriguing 
finding  of  this  study  is  that  overexpressing  of  an  inactive  form  of  ERK1  or  ERK2  alone  is 
sufficient  to  block  the  actions  of  endogenous  ERKs  in  response  to  flow  shear  stress  in  bone  cell 
proliferation. 

Specific  Objective  2:  Evaluate  changes  in  gene  expression  of  those  genes  known  to 
be  involved  in  mediating  loading  signal  using  in-house  microarray. 

In  order  to  study  the  genes  involved  in  the  mechanical  signaling  pathway,  we  applied  the 
microarray  technique.  Our  hypothesis  is  that  the  expression  level  of  the  genes  involved  in  the 
mechanical  signaling  would  be  higher  in  the  B6  mice  as  compared  to  C3H  mice  since  B6  mice 
are  responsive  to  mechanical  stress  whereas  C3H  mice  do  not  respond  to  mechanical  loading. 
Thus,  to  study  the  genetic  basis  for  difference  in  loading  signal  in  response  to  shear  stress 
between  C3H  and  B6  mice  strains,  we  identified  the  genes  in  which  the  expression  level 
increased  more  than  two-fold  when  subjected  to  shear  stress  as  compared  to  no  stress  in  these  two 
strains.  The  bone  cells  isolated  from  B6  and  C3H  mice  were  subjected  to  fluid  flow  shear  stress 
of  20  dynes/cm2  for  30  minutes.  After  4  hours,  RNA  was  extracted  and  microarray  was 
performed. 

Significantly  Increased  or  Decreased  Gene  Expression  from  B6  Osteoblasts  After 
Sheer  Stress  (Tables  1  &  2) 

Table  1.  List  of  genes  in  which  expression  level  increased  more  than  two-fold  after 
application  of  shear  stress  as  compared  to  control  cells  isolated  from  B6  mice.  The  ratios 
for  the  microarray  were  determined  by  using  Genespring  software. 


Accession  # 

Fold 

Change 

C81354 

6.021867 

XI 5848 

3.62071 

J04113 

3.536071 

AF013170 

3.508688 

TCI 59688 

3.315779 

X05010 

3.117075 

AU024669 

3.105589 

M69293 

3.028576 

E04743 

2.992537 

X57796 

2.960955 

J00370 

2.927765 

J05265 

2.894687 

U64331 

2.887784 

AF00901 1 

2.781368 

AB015978 

2.778422 

Gene  Name 

Retinoic  acid  receptor  gamma 

Mouse  thyroid  hormone  receptor  (N UR/77) 

TNF  related  ligand  TRANCE 

Mus  musculus  IGFBP5  mRNA,  complete  cds 

Colony  stimulating  factor  alpha  (TNFa)  gene 

Mouse  Id2  protein  (Id-2) 

Mouse  IL-1  alpha 

Mouse  55-kda  tumor  necrosis  factor  receptor 
Mouse  c-fos  gene;  cellular  homolog  to  viral  oncogene 
Interferon  gamma  receptor 
Mouse  Osteoprotegerin  (OPG) 
axin 

Mouse  Oncostatin  M  receptor  beta 


U81451 

M28021 

C78332 

M59929 

U58503 

E01057 

D49438 

X56848 

AF053713 

AJ009862 

L35303 

L0253F1 1-3 

L28108 

L0065H09-3 

TC201479 

TCI  75431 

D89628.1 

BB500047 

U39060 

U36384 

TCI  57862 

U85259 

XI 5202 

M89798 

U51001 

SU78076 

SJ00357 

U78048 

D1 6250.1 

AF019048 

TC200632 

AF068615 

NM_0 10703 

X76401 

J04069 

SU96809 

TCI 77706 

C78942 

C80859 

V00727 

D17630 

TCI  87351 

D31942 

AF1 79369 

C85074 

TC202099 

E03515 

J04766.1 


2.773452  Estrogen  receptor  beta  (ESTRB) 

2.752648  Mouse  Homeo  domain  protein  (HOX  1.3) 

2.727276 

2.726371  Signal-trans.  guanine  nuleotide-binding  protein  (GNA01) 

2.682652  Keratinocyte  growth  factor/fibroblast  growth  factor-7 
2.672017  Mouse  interleukin-1  (IL-1 )  precursor 
2.64229  Mouse  25-hydroxyvitamin  D3  24-hydroxylase 
2.610531  Bone  morphogenetic  protein  4  (BMP-4) 

2.600324  Mouse  Osteoprotegerin  ligand 
2.584103  Transforming  growth  factor-beta  1 
2.571284  TNF  receptor  associated  factor  2  (TRAF2) 

2.566142  Mus  musculus  orphan  nuclear  hormone  receptor  (CAR)  gene,  complete  sequence 
2.553446  Mouse  PTH/PTHrP 

2.529236  Homo  sapiens  ribosomal  protein  S21  (RPS21),  mRNA 

2.526702  Homo  sapiens  MEBP-1  mRNA  for  MAPK-ERK  binding  protein-1,  complete  cds 

2.519852 

2.506003  Vascular  endothelial  growth  factor  D  (VEGF-D) 

2.498375  beta-catenin 

2.49294  Glucocorticoid  receptor  interacting  protein  1  (GRIP1) 

2.488445  Mouse  twist-related  bHLH  protein  Dermo-1 

2.438857  Homo  sapiens  mRNA;  cDNA  DKFZp586C1620  (from  clone  DKFZp586C1620);  partial  cds 
2.437383  Estrogen  related  receptor  alpha  (ESTRRA) 

2.428389  Fibronectin  Receptor  beta-chain  (VLA5-homolog.) 

2.414296  Mouse  WNT-5a 

2.398422  DLX-1  gene 

2.394921  Mouse  sepiaterin  reductase  gene 

2.394785  Mouse  alpha  amylase-2  gene 

2.38784  Bone  morphogenetic  protein  type  II  receptor  BRK-3 

2.37863  Mouse  BMP  receptor 

2.370698  Mouse  RANKL 

2.360679  Homo  sapiens  CCR4-NOT  transcription  complex,  subunit  8  (CNOT8),  mRNA 
2.358135  Ciliary  neutotrophic  factor  receptor  alpha  precursor 
2.355529  Lef  1 

2.3401 1 1  Tumor  necrosis  factor  receptor  2 
2.338559  Insulin-like  growth  factor  II  (IGF-II) 

2.333815  Mouse  chromatin  structural  protein  homolog  gene 

2.332596 

2.326618 

2.321968 

2.315034  c-fos  oncogene 

2.31 364  Mouse  mRNA  for  interleukin-8  receptor 

2.313193  Homo  sapiens  actin  related  protein  2/3  complex,  subunit  4  (20  kD)  (ARPC4),  mRNA 
2.295312  Mouse  Oncostatin  M 

2.28034  Insuline-like  growth  factor  binding  protein  5  protease 

2.27293 

2.258783 

2.257122  Mouse  interleukin  6  receptor  protein 
2.254782  Plasminogen 


L0029F01-3 

2.248046 

Y15800 

2.242617 

L0051E10-3 

2.22982 

C76615 

2.224257 

TCI 88980 

2.223752 

TCI 86487 

2.210573 

TC161819 

2.210151 

LI 2447 

2.194287 

TCI 86750 

2.191918 

M92416 

2.191191 

TC1661 14 

2.186043 

M29464 

2.179255 

TC160957 

2.176756 

TC162199 

2.171346 

TCI 86822 

2.155771 

C81425 

2.154317 

X53802 

2.150552 

TC176196 

2.149525 

SM95800 

2.14905 

TC160277 

2.138307 

TC201280 

2.131578 

J04953 

2.129044 

TCI  86937 

2.125651 

AF067191 

2.124936 

TCI  8801 8 

2.116157 

TCI  74028 

2.107252 

TC160410 

2.089787 

TCI  75783 

2.08736 

H3147A10 

2.081017 

TCI  87978 

2.074785 

AF1 20489 

2.073418 

M63650 

2.071397 

TCI  74003 

2.070687 

TCI  74701 

2.069076 

TC200417 

2.060727 

TCI 73360 

2.057717 

AF075717 

2.044671 

TCI  99965 

2.043639 

TC162168 

2.041731 

C80763 

2.038737 

AU022524 

2.028801 

X58636.1 

2.025792 

H3120A04 

2.021601 

C86855 

2.018469 

L32751 

2.013074 

AA1 19293 

2.012701 

TCI  86472 

2.00335 

TC201234 

2.000614 

Mouse  mRNA  for  G-protein  coupled  receptor  kinase  6-D 
Mus  musculus  fibrillarin  (Fbl),  mRNA 


Mus  musculus  (clone  Clebp-1)  high  mobility  group  1  protein  (HMG-l)mRNA, complete  cds 

Insuline-like  growth  factor  binding  protein  5  (IGFBP5) 

Mus  musculus  protein  tyrosine  phosphatase  4a1  (Ptp4a1),  mRNA 
Fibroblast  growth  factor  (FGF6) 

Mouse  mRNA  for  ELP3,  complete  cds 

Mouse  platelet-derived  growth  factor  a  chain  (PDGA) 

Mus  musculus  PL-6  (PI6),  mRNA 

Homo  sapiens  splicing  factor  (CC1.3)  (CC1.3),  mRNA 

Mouse  mRNA  for  interleukin-6  receptor 

Mouse  Myogenin  gene 

Mus  musculus  15  kDa  selenoprotein  (Sep15)  mRNA,  complete  cds 
Mus  musculus  D-dopachrome  tautomerase  gene,  complete  cds 
Gelsolin  gene 

Mus  musculus  DBA/2J  delta  proteasome  subunit  gene,  complete  cds 
Fibroblast  growth  factor  1  .a  gene,  exon  1 

Mus  musculus  spermidine/spermine  Nl-acetyl  transferase  (Sat),  mRNA 
Mus  musculus  serum  and  glucocorticoid-dependent  protein  kinase  (Sgk)mRNA 
Mus  musculus  hyaluronidase  2  (Hyal2),  mRNA 

Mus  musculus  recombining  binding  protein  suppressor  of  hairless(Drosophila)(Rbpsuh) 
Growth  hormone  receptor/growth  hormone-binding  protein 
Mouse  M-twist  gene 

Mus  musculus  kidney  predominant  protein  NCU-G1  (NCU-G1),  mRNA 
Homo  sapiens  JTV1  gene  (JTV1),  mRNA 

Human  gene  fordihydrolipoamide  succinyltransferase,  complete  cds  (exon  1-15) 

Homo  sapiens  actin  related  protein  2/3  complex,  subunit  5  (16  kD)  (ARPC5),  mRNA 
Transforming  growth  factor  beta  1 -induced  factor  2 

Homo  sapiens  cDNA  FLJ14219  fis,  clone  NT2RP3003800,  highly  similar  to  Rattus  norveg 
Homo  sapiens  CGI-147  protein  (LOC51651),  mRNA 


(LEF1)  lymphoid  enhancer  binding  factor  1 

Mus  musculus  scmhl  mRNA  for  sex  comb  on  midleg  homolog  protein,  complete  cds 

Mouse  (Clone  Ml)  GTPase  (RAN) 

small  inducible  cytokine  A6 

Mouse  t-complex  protein  (Tcp-1  x)  mRNA,  3’  end 

Rattus  norvegicus  lamina  associated  polypeptide  1C  (LAP1C)  mRNA,  complete  cds 


Table  2.  List  of  genes  in  which  expression  level  decreased  more  than  two-fold  after 
application  of  shear  stress  as  compared  to  control  cells  isolated  from  B6  mice.  The  ratios 
for  the  microarray  were  determined  by  using  Genespring  software. 


Accession  # 

TC163213 
TCI  86806 
C79115 
C87796 
TCI  77345 
TCI  931 80 
TC167104 
C80790 
C79008 
TCI  60201 
TCI  76539 
TCI  90432 
TCI  73689 
TCI  76094 
TCI  59792 
L0298F03-3 
AU019373 
C86899 
TC202698 
TCI 58511 
L0069G07-3 
TCI 74040 
TCI  88831 
TCI  78609 
TC204210 
TC208948 
AU024699 
C81530 
AU022331 
TC204243 
TCI 75962 
TCI  89875 
TCI 62398 
L0223G04-3 
TC161135 
TCI 63534 
TCI 90975 
AU020234 
C87530 
L0226E04-3 


Fold 

Change  Gene  Name 

6.323833 

4.749427  Homo  sapiens  apoptosis-related  protein  PNAS-3  (PNAS-3)  mRNA,  partial  cds 

4.397889 

4.207526 

4.145114  Homo  sapiens  hypothetical  protein  FLJ1 0761  (FLJ1 0761),  mRNA 

3.797844  Mus  musculus  tyrosine  hydroxylase  (Th),  mRNA 

3.656022 

3.562263 

3.479888 

3.372998  Homo  sapiens  mRNA;  cDNA  DKFZp586A0722  (from  clone  DKFZp586A0722) 

3.316524  Mus  musculus  small  GTP-binding  protein  RAB25  (Rab25)  gene,  complete  cds 
3.315813 

3.286380  Mus  musculus  putative  CCAAT  binding  factor  1  (mCBF)  mRNA,  mCBFI ,  complete  cds 

3.229843  Homo  sapiens  KIAA0205  gene  product  (KIAA0205),  mRNA 

3.224493  Homo  sapiens  ribosomal  protein  S9  (RPS9),  mRNA 

3.224286 

3.178850 

3.094014 

3.087886  Homo  sapiens  zinc  finger  protein  198  (ZNF198),  mRNA 

3.063481  Homo  sapiens  mRNA;cDNA  DKFZp586P0123  (from  clone  DKFZp586P0123);partial  cds 

3.044982  Homo  sapiens  hypothetical  protein  FLJ10579  (FLJ10579),  mRNA 

3.028745 

3.027351  Homo  sapiens  small  nuclear  ribonucleoprotein  polypeptide  F  (SNRPF),  mRNA 

3.017881  Mus  musculus  nuclear  RNA  export  factor  1  homolog  (S.  cerevisiae)  (Nxfl),  mRNA 

3.013675  Homo  sapiens  B-cell  CLL/lymphoma  9  (BCL9),  mRNA 

3.009861  Mus  musculus  Janus  kinase  2  (Jak2),  mRNA 

2.999606 

2.965025 

2.935256 

2.907804  Homo  sapiens  mRNA  for  KIAA1 524  protein,  partial  cds 

2.872361  Homo  sapiens  clone  CTB-10G5,  complete  sequence 

2.871053  Homo  sapiens  hypothetical  protein  FLJ11110  (FLJ11110),  mRNA 

2.857831  Homo  sapiens  insulin  receptor  tyrosine  kinase  substrate  (LOC55971),  mRNA 

2.835623 

2.828019  Homo  sapiens  cDNA  FLJ 14228  fis,  clone  NT2RP30041 48 
2.797985 

2.791493  Homo  sapiens  centrosomal  P4.1 -associated  protein  (CPAP)  mRNA,  complete  cds 
2.767904  Homo  sapiensTRPM-2, cytosolic  epoxide  hydrolase, nicotinic  acetylcholine  receptor  alph; 
2.755052 

2.734502  Homo  sapiens  LIM  domain  only  7  (LM07),  mRNA 


TC202077 
TCI  91 546 
C80806 
TCI  78272 
C81384 
L0051A02-3 
H3129D12 
TCI  63806 
TCI 5721 9 
C81435 
TCI 871 87 
TC189163 
C78157 
L0290G10-3 
TCI 57805 
TCI  73604 
TC2 10779 
TCI  59887 
TCI  88500 
TCI 60230 
TCI  74370 
NM-008513 
TCI 78738 
TCI 76946 
TCI 89678 
TC161711 
TCI  86798 
TCI  88265 
TC206463 
C86942 
TC204786 
TCI  741 82 
S67455 
TC163121 
TCI 761 86 
L0291D11-3 
L0013H08-3 
TCI 90732 
TCI  92702 
EST03087 
TCI  77370 
TCI  76521 
TCI  59038 
TC201917 
TC191320 
AU 0224 14 
TC160509 
C85917 


2.674004 

2.673286  Homo  sapiens  cDNA  FLJ12145  fis,  clone  MAMMA1 000395 
2.666355 

2.665531  Homo  sapiens  cDNA  FLJ1 3303  fis, clone  OVARC1 001 372, highly  similar  to  Homo  sapien 
2.663685  Homo  sapiens  KIAA0042  gene  product  (KIAA0042),  mRNA 
2.622076  Homo  sapiens  heterogeneous  nuclear  ribonucleoprotein  A0  (HNRPA0),  mRNA 
2.609591 

2.601031  Homo  sapiens  heat  shock  transcription  factor  2  binding  protein  (HSF2BP),  mRNA 
2.574551  Homo  sapiens  eukaryotic  protein  synthesis  initiation  factor  mRNA,  complete  cds 
2.564388 

2.547066  Homo  sapiens  mRNA;  cDNA  DKFZp761M2324  (from  clone  DKFZp761M2324) 

2.544920  Homo  sapiens  cDNA:  FLJ23321  fis,  clone  HEP12396 

2.535713 

2.516292  Homo  sapiens  KIAA0738  gene  product  (KIAA0738),  mRNA 
2.512923  Homo  sapiens  GTP-binding  protein  (NGB),  mRNA 

2.511347  Homo  sapiens  loss  of  heterozygosity,  1 1 ,  chromosomal  region  2,  gene  A  (LOH11CR2A) 
2.507478  Homo  sapiens  KIAA0513  gene  product  (KIAA0513),  mRNA 
2.492242  Homo  sapiens  KIAA0005  gene  product  (KIAA0005),  mRNA 
2.489737  Homo  sapiens  mRNA  for  repressor  protein,  partial  cds 

2.480316  Homo  sapiens  mRNA;  cDNA  DKFZp566J151(from  clone  DKFZp566J151);complete  cds 
2.479888  Mus  musculus  solute  carrier  family  34  (sodium  phosphate),  member  2  (Slc34a2),  mRNA 
2.473082  LRP5  #6 
2.470564 

2.458136  Homo  sapiens  hypothetical  protein  FLJ20241  (FLJ20241),  mRNA 

2.452885  Homo  sapiens  mRNA  for  KIAA0693  protein,  partial  cds 

2.422831  Homo  sapiens  nesca  protein  (NESCA),  mRNA 

2.416530  Homo  sapiens  hypothetical  protein  FLJ10488  (FLJ10488),  mRNA 

2.399167  Homo  sapiens  Sec31  protein  mRNA,  complete  cds 

2.385723 

2.350195 

2.340445  Homo  sapiens  SWI/SNFrelated, matrix  associated, actin  dependent  regulator  of  chromatii 
2.333507  Mus  musculus  platelet-activating  factor  acetylhydrolase,  isoform  1b,(Pafah1b3),  mRNA 
2.324881  Osteocalsin 

2.321306  Homo  sapiens  cDNA  FLJ12969  fis,  clone  NT2RP2005841 
2.298626  Homo  sapiens  LIM  and  cysteine-rich  domains  1  (LMCD1),  mRNA 
2.279614 

2.278326  Homo  sapiens  eukaryotic  translation  elongation  factor  1  gamma  (EEF1G),  mRNA 

2.274502  Homo  sapiens  cDNA:  FLJ21020  fis,  clone  CAE06067 

2.261900 

2.260475 

2.242102  Human  Tcr-C-delta  gene, exons  1-4;  Tcr-V-delta  gene, exons  1-2  Tcr-C-alpha  gene,  exor 

2.238748  Homo  sapiens  KIAA0615  gene  product  (KIAA0615),  mRNA 

2.221344  Homo  sapiens  mRNA;  cDNA  DKFZp761E212  (from  clone  DKFZp761E212) 

2.219296  Homo  sapiens  cDNA  FLJ1 2783  fis,  clone  NT2RP2001 876 

2.217481  Rattus  norvegicus  gene  forT  cell  receptor  eta  chain,  partial  cds 

2.214434  Homo  sapiens  KIAA0940  protein  (KIAA0940),  mRNA 

2.214421  Homo  sapiens  diaphorase  (NADH)  (cytochrome  b-5  reductase)  (DIA1),  mRNA 

2.205803  Homo  sapiens  spinocerebellar  ataxia  1 


TCI  90467 
TC203547 
TC160176 
TCI 61 597 
TC202599 
TCI 90059 
TCI  57266 
TCI  99932 
TC161470 
TCI 63446 
TCI  78608 
L0261D05-3 
TCI 90403 
TCI 67570 
TCI  57243 
TCI  58841 
TCI 60307 
TCI  76029 
L0221A07-3 
C81146 
TCI  60464 
TC202524 
TCI  89397 
TCI  89007 
TCI  88371 
TCI  88050 
C78068 
TCI  63941 
TC204346 
AU019250 
L0215A02-3 
TC200441 
TCI  91 257 
TCI  87210 
TCI 7381 9 
C76739 
AU019952 
AU041136 
TCI  88786 
AU024767 
L0290B05-3 
TCI 60773 
AUO 19202 
TC204088 
AU022194 
TC163812 
TCI 64591 
TC176191 


2.199626  Homo  sapiens  cDNA  FLJ 13645  fis,  clone  PLACE1 01 1310 
2.198232  Homo  sapiens  cDNA  FLJ  12642  fis,  clone  NT2RM4001 965 
2.198150  Homo  sapiens  mRNA  for  KIAA1470  protein,  partial  cds 
2.194239 

2.181680  Homo  sapiens  KIAA0244  protein  (KIAA0244),  mRNA 

2.174398  Homo  sapiens  transmembrane  4  superfamily  member(tetraspan  NET-7)  (NET-7), mRNA 
2.169345  Mus  musculus  macrophage  galactose  N-acetyl-galactosamine  specific  lectin(Mgl),mRN/ 
2.167705  Homo  sapiens  MADS  box  transcription  enhancer  factor2,  polypeptide  B  (MEF2B),mRN4 
2.163998  Homo  sapiens  mRNA;  cDNA  DKFZp434E0121  (from  clone  DKFZp434E0121) 

2.158835  Mus  musculus  murinoglobulin  1  (Mugl),  mRNA 

2.158464  Homo  sapiens  8q22.1  region  and  MTG8  (CBFA2T1)  gene,  partial  cds 

2.158177  Homo  sapiens  similar  to  phosphorylase  kinase,  alpha  2  (liver)  (H.  sapiens)  (LOC63600) 

2.152387  Homo  sapiens  G-protein-coupled  receptor  induced  protein  GIG2  (GIG2)  mRNA 

2.146512  Mus  musculus  citrin  (Slc25a1 3)  mRNA,  complete  cds 

2.139440  Homo  sapiens  myosin  1C  (MYOIC),  mRNA 

2.133129  Homo  sapiens  cDNA  FLJ  13924  fis,  clone  Y79AA1 000540 

2.131873  Homo  sapiens  isocitrate  dehydrogenase  3  (NAD+)  alpha  (IDH3A),  mRNA 

2.128533  Homo  sapiens  mRNA;  cDNA  DKFZp434L1850(from  clone  DKFZp434L1 850); partial  cds 

2.125280 

2.119302 

2.116825  Homo  sapiens  mRNA-associated  protein  mrnp41  mRNA,  complete  cds 

2.112259  Human  TFIID  subunits  TAF20  and  TAF15  mRNA,  complete  cds 

2.096525  Homo  sapiens  hypothetical  protein  (FLJ20323),  mRNA 

2.094517  Mus  musculus  RAD50  homolog  (S.  cerevisiae)  (Rad50),  mRNA 

2.094454  Human  clone  23721  mRNA  sequence 

2.093248  Homo  sapiens  MDS017  (MDS017)  mRNA,  complete  cds 

2.078961 

2.074859 

2.066088  Homo  sapiens  cDNA:  FLJ22573  fis,  clone  HSI02387 
2.063553  Homo  sapiens  T-box  19  (TBX19),  mRNA 
2.058678 

2.047412  Homo  sapiens  guanine  nucleotide  binding  protein  (G  protein),  (GNAI3),  mRNA 

2.044590  Human  IGF-I  mRNA  for  insulin-like  growth  factor  I 

2.040991  Homo  sapiens  cDNA  FLJ  13872  fis,  clone  THYR01 001322 

2.035938  Homo  sapiens  thyroid  receptor  interactor  (TRIP3)  mRNA,  3'  end  of  cds 

2.032122  Mus  musculus  macrophage  C-type  lectin  (Mpcl),  mRNA 

2.032041  Homo  sapiens  mRNA  for  KIAA0597  protein,  partial  cds 

2.027832  Bacteriophage  lambda,  complete  genome 

2.026039  Rattus  norvegicus  NADH/NADPH  mitogenic  oxidase  subunit  p65-mox  mRNA 

2.024597  Mus  musculus  mRNA  for  PC3B  protein 

2.023050 

2.020360  M. musculus  DNA  3'flanking  minisatellite  transgene  HOC 
2.015859  Homo  sapiens  MSTP046  mRNA,  complete  cds 
2.015016  Homo  sapiens  HYA22  protein  (HYA22),  mRNA 
2.010488 

2.010332  Homo  sapiens  cDNA  FLJ  10366  fis,  clone  NT2RM2001420 
2.009105  M. musculus  mRNA  for  intestinal  tyrosine  kinase 
2.008425  Homo  sapiens  BM022  protein  (BM022),  mRNA 


AU 024701  2.008356 

TC203734  2.000904  Homo  sapiens  mRNA  for  KIAA1 302  protein,  partial  cds 

C78087  2.000170 

As  shown  in  Table  1,  there  are  109  genes  that  show  more  than  two-fold  increase  in 
expression  in  cells  subjected  to  fluid  flow  shear  stress  as  compared  to  the  non-stressed  cells  in  B6 
mice.  Table  2  shows  that  there  are  139  genes  that  show  more  than  two-fold  decrease  in 
expression  in  stressed  cells  as  compared  to  the  non-stressed  cells  in  B6  mice. 

Significantly  Increased  or  Decreased  Gene  Expression  in  C3  Osteoblasts  After  Sheer 
Stress  (Tables  3  &  4) 

Table  3.  List  of  genes  in  which  expression  level  increased  more  than  two-fold  after 
application  of  shear  stress  as  compared  to  control  cells  isolated  from  C3H  mice.  The 
ratios  for  the  microarray  were  determined  by  using  Genespring  software. 


Accession  # 

Fold 

Change 

M69293 

5.342443 

AF1 79369 

4.480995 

M63650 

4.314368 

D16250.1 

4.153827 

M89798 

3.86404 

U85259 

3.828085 

AF056187 

3.623736 

AF126159 

3.500622 

X57413 

3.49456 

L27424 

3.46014 

U39060 

3.454962 

U58503 

3.376605 

AJ009862 

3.354771 

D 17630 

3.210925 

U67610 

3.206105 

X05010 

3.13227 

AB015978 

3.128478 

M92416 

3.113992 

AB009993 

3.08771 

M97017 

3.059502 

L25602.1 

3.019166 

X67348 

3.012249 

XI 5202 

2.947969 

J00370 

2.937989 

TCI 77706 

2.931965 

U39545 

2.926427 

M59929 

2.915431 

U70429 

2.847221 

U81451 

2.841179 

Gene  Name 

Mouse  Id2  protein  (Id-2) 

Insuline-like  growth  factor  binding  protein  5  protease 
Mouse  M-twist  gene 
Mouse  BMP  receptor 
Mouse  WNT-5a 

Estrogen  related  receptor  alpha  (ESTRRA) 

Insulin-like  growth  factor  I  receptor 
Mouse  big  MAP  kinase  la  (BMK1),  mRNA 
Transforming  growth  factor-beta  2 
Mouse  metalloproteinase  inhibitor  TIMP-3) 
Glucocorticoid  receptor  interacting  protein  1  (GRIP1) 
Keratinocyte  growth  factor/fibroblast  growth  factor-7 
Transforming  growth  factor-beta  1 
Mouse  mRNA  for  interleukin-8  receptor 
Fibroblast  growth  factor  1  (FGF-1) 

Colony  stimulating  factor  alpha  (TNFa)  gene 
Mouse  Oncostatin  M  receptor  beta 
Fibroblast  growth  factor  (FGF6) 

Collagen  A1(V) 

Osteogenic  protein-2  (OP-2) 

Mouse  morphogenetic  protein  2  (BMP-2)  gene 
COL10A1  gene  for  alpha  -  collagen  type  x 
Fibronectin  Receptor  beta-chain  (VLA5-homolog.) 
Mouse  c-fos  gene;  cellular  homolog  to  viral  oncogene 

Bone  morphogenetic  protein  8b  (BMP8b) 

Signal-trans.  guanine  nuleotide-binding  protein  (GNA01) 
lnterleukin-4  induced  gene-1  (FIG1) 

Estrogen  receptor  beta  (ESTRB) 


J04953 

2.811038 

AB006034 

2.796519 

L32751 

2.788244 

M29464 

2.778772 

L35303 

2.759914 

AF020681 

2.758174 

TCI 75431 

2.745054 

D89628.1 

2.738802 

XI 4759 

2.72283 

U64331 

2.718682 

X57796 

2.678408 

D63644 

2.6558 

E04743 

2.655764 

TC200162 

2.654953 

XI 5848 

2.650977 

M28021 

2.638262 

TC173131 

2.605621 

AB021 228.1 

2.602641 

TC173318 

2.540124 

NM_0 10703 

2.507729 

J04113 

2.482183 

TCI 88890 

2.457797 

TCI  87062 

2.448608 

Perkin  Elmer 

2.425088 

C86855 

2.401036 

SU78076 

2.395678 

J05265 

2.391902 

U36384 

2.39182 

SM95800 

2.381928 

X62622 

2.374983 

TCI 58549 

2.365563 

TC161819 

2.363219 

TCI 90585 

2.352407 

J 04069 

2.350386 

L0227G08-3 

2.347008 

LI  5436 

2.336559 

TC160410 

2.316149 

TCI 86487 

2.30955 

H3120A04 

2.303708 

AF067191 

2.290257 

L28108 

2.288727 

TCI 87351 

2.267528 

L0011C07-3 

2.266767 

TC201280 

2.252388 

AA253928 

2.252074 

TC160692 

2.249222 

H3124G02 

2.242975 

C78676 

2.236116 

Gelsolin  gene 

Mouse  25-hydroxyvitamin  D3  1  alpha-hydroxylase 
Mouse  (Clone  Ml)  GTPase  (RAN) 

Mouse  platelet-derived  growth  factor  a  chain  (PDGA) 

TNF  receptor  associated  factor  2  (TRAF2) 

Mouse  core  binding  factor  alphal  sub  unit  isoform  (Cbfa-1) 

Vascular  endothelial  growth  factor  D  (VEGF-D) 

Mouse  Homeo  box  (HOX-7.1) 

Mouse  Osteoprotegerin  (OPG) 

Mouse  55-kda  tumor  necrosis  factor  receptor 
ARNT2  (F2#335) 

Mouse  IL-1  alpha 

Homo  sapiens  ribosomal  protein  L39  (RPL39),  mRNA 

Retinoic  acid  receptor  gamma 

Mouse  Homeo  domain  protein  (HOX  1.3) 

Homo  sapiens  eukaryotic  translation  initiation  factor  2, subunit  2(beta,  38kD  )(EIF2S2) 

Membrane-type  3  matrix  metalloproteinase 

Homo  sapiens  transmembrane  trafficking  protein  (TMP21),  mRNA 

Lef  1 

Mouse  thyroid  hormone  receptor  (N UR/77) 

Mus  musculus  calponin  1  (Cnnl),  mRNA 
Homo  sapiens  6.2  kd  protein  (LOC54543),  mRNA 
EST  91  from  A.Thaliana 

Mus  musculus  scmhl  mRNA  for  sex  comb  on  midleg  homolog  protein,  complete  cds 

Mouse  sepiaterin  reductase  gene 

Interferon  gamma  receptor 

Mouse  twist-related  bHLH  protein  Dermo-1 

Mouse  Myogenin  gene 

(TIMP-2)  Tissue  inhibitor  of  metalloproteinases 


Insulin-like  growth  factor  II  (IGF-II) 

Homo  sapiens  clone  82F9,  complete  sequence 
Isoform  of  TGF-b  type  II  receptor 

Mus  musculus  serum  and  glucocorticoid-dependent  protein  kinase  (Sgk)  mRNA 
Mus  musculus  (clone  Clebp-1)  high  mobility  group  1  protein  (HMG-1)  mRNA 

Fibroblast  growth  factor  la  gene,  exon  1 
Mouse  PTH/PTHrP 

Homo  sapiens  actin  related  protein  2/3  complex,  subunit  4  (20  kD)  (ARPC4),  mRNA 

Mus  musculus  D-dopachrome  tautomerase  gene,  complete  cds 
endothelial  monocyte 

Homo  sapiens  myeloid  leukemia  factor  2  (MLF2),  mRNA 


AF026305 

2.230761 

Y15800 

2.227902 

TCI 60740 

2.21674 

C78942 

2.21332 

C87911 

2.21322 

X81582 

2.210164 

TCI  91 001 

2.210023 

AF 126063 

2.209931 

C81495 

2.208941 

AA1 19293 

2.203506 

TC159712 

2.199665 

TCI  59899 

2.198992 

D31942 

2.197945 

TCI  73358 

2.197787 

C81364 

2.179308 

TCI  60481 

2.174184 

TCI  59868 

2.173848 

TCI  86472 

2.169085 

TCI  86937 

2.168598 

AU022429 

2.168419 

C81299 

2.167638 

TC161269 

2.166696 

C76825 

2.154873 

TCI  90469 

2.149796 

TC201640 

2.149553 

AU024669 

2.149365 

L0253F11-3 

2.146557 

L0065E08-3 

2.146445 

C78257 

2.145905 

AF00901 1 

2.145423 

C87927 

2.139591 

TCI  74006 

2.135762 

AF013170 

2.135479 

X56848 

2.13524 

AA067193 

2.131849 

TC174140 

2.128732 

TCI 60359 

2.125901 

D49438 

2.125615 

TC200023 

2.119359 

TCI  74370 

2.117575 

C86037 

2.107628 

BB500047 

2.099311 

L0041A01-3 

2.090664 

C78966 

2.077118 

L0065H09-3 

2.074872 

C78958 

2.068447 

TC189109 

2.065315 

TC181404 

2.061944 

Mouse  Zinc  finger  transcription  factor  GLI 

Mouse  mRNA  for  G-protein  coupled  receptor  kinase  6-D 

Homo  sapiens  translocation  protein  1  (TLOC1),  mRNA 

Mus  musculus  SIL,  MAP_17,  CYP_a,  SCL  &  CYP_b  genes 
Insulin-like  growth  factor  binding  protein-4 
Mus  musculus  transcription  factor  LRG-21  mRNA,  complete  cds 
Connective  tissue  growth  factor-like  protein  precursor 

small  inducible  cytokine  A6 

Mus  musculus  RNA  polymerase  1-1  (40  kDa  subunit)  (Rpol-1),  mRNA 
Mouse  Oncostatin  M 

Homosapienssignal  sequence  receptor, beta(translocon-associated  protein  beta)(SSR2) 

Mus  musculus  calcium-binding  protein  Cab45a  mRNA,  complete  cds 

Mus  musculus  necdin  (Ndn),  mRNA 

Mouse  t-complex  protein  (Tcp-1  x)  mRNA,  3'  end 

Mus  musculus  DBA/2J  delta  proteasome  subunit  gene,  complete  cds 


Homo  sapiens  hypothetical  protein  FLJ20272  (FLJ20272),  mRNA 
Mus  musculus  transient  receptor  potential-related  protein  (ChaK),  mRNA 

Mus  musculus  arginine  methyltransferase  (Prmt2)  mRNA,  complete  cds 

Mus  musculus  orphan  nuclear  hormone  receptor  (CAR)  gene,  complete  sequence 

Human  mRNA  for  NADP  dependent  leukotriene  b4  12-hydroxydehydrogenase, partial  cds 
axin 

Homo  sapiens  asparaginyl-tRNA  synthetase  (NARS),  mRNA 

TNF  related  ligand  TRANCE 

Bone  morphogenetic  protein  4  (BMP-4) 

UDP-glucose  dehydrogenase 

Mus  musculus  Cdc42  GTPase-inhibiting  protein  (Cdgip-pending),  mRNA 
Homo  sapiens  clone  RP1 1-486122,  complete  sequence 
Mouse  25-hydroxyvitamin  D3  24-hydroxylase 
Mus  musculus  ribosomal  protein  L27  (Rpl27),  mRNA 

Mus  musculus  solute  carrier  family  34  (sodium  phosphate),  member  2  (Slc34a2),  mRNA 
beta-catenin 


Homo  sapiens  ribosomal  protein  S21  (RPS21),  mRNA 

Mus  musculus  transcription  elongation  factor  B  (Sill),  polypeptide  3  (1 10kD)  (Tceb3) 


TCI 60525 

2.06007 

L0029F01-3 

2.059959 

L0008A03-3 

2.048002 

TC201527 

2.047403 

TC187519 

2.047178 

TCI 57688 

2.046966 

TCI 73042 

2.03422 

U78048 

2.028893 

AU040912 

2.02622 

TCI 57249 

2.020829 

TCI  60904 

2.015319 

C76812 

2.013664 

TCI  60230 

2.011371 

SJ00357 

2.00887 

TC159122 

2.007755 

C81549 

2.005677 

TC157612 

2.002072 

X76401 

2.001228 

Homo  sapiens  cDNA:  FLJ21894  fis,  clone  HEP03434 


Homo  sapiens  from  HeLa  cyclin-dependent  kinase  2  interacting  protein  (CINP),  mRNA 

Mus  musculus  hematological  and  neurological  expressed  sequence  1  (Hnl),  mRNA 
Homo  sapiens  NADH-ubiquinone  dehydrogenase  1  beta  subcomplex  mRNA,  complete  cc 
Bone  morphogenetic  protein  type  II  receptor  BRK-3 

Homo  sapiens  mRNA  for  KIAA0622  protein,  partial  cds 


Homo  sapiens  mRNA;  cDNA  DKFZp566J151  (from  clone  DKFZp566J151);  complete  cds 
Mouse  alpha  amylase-2  gene 

Homo  sapiens  clone  RP1 1-359J14,  complete  sequence 

Mus  musculus  thioredoxin  reductase  1  (Txnrdl),  mRNA 
Tumor  necrosis  factor  receptor  2 


Table  4.  List  of  genes  in  which  expression  level  decreased  more  than  two-fold  after 
application  of  shear  stress  as  compared  to  control  cells  isolated  from  C3H  mice.  The 
ratios  for  the  microarray  were  determined  by  using  Genespring  software. 


Accession  # 

Fold 

Change 

L0003B11-3 

7.251633 

C81354 

5.150292 

AU024549 

4.566669 

C86942 

4.135808 

TCI 61 706 

3.814792 

AU023208 

3.674891 

TC203078 

3.630153 

TCI  69364 

3.333058 

AU019631 

3.262583 

TCI  58443 

3.253914 

AU040173 

3.189639 

TCI 73754 

3.153296 

TCI  87226 

3.104026 

AU040981 

2.964905 

TC171051 

2.935193 

C78087 

2.929459 

EST03087 

2.875889 

TCI  62350 

2.870055 

TC160175 

2.781838 

Gene  Name 

Homosapiens  FRA3B  common  fragile  region, diadenosine  triphosphate  hydrolase  (FHf 


Mus  musculus  guanine  nucleotide  binding  protein,  alpha  14  (Gna14),  mRNA 


Mus  musculus  DAZ-like  putative  RNA  binding  protein  mRNA,  complete  cds 
Mus  musculus  EST  from  clone  1498755,  3'  end 


Mus  musculus  UDP-N-acetyl-alpha-D-galactosamine:polypeptideN-  3  (Galnt3),  mRNA 
Mus  musculus  serine/threonine  kinase  receptor  associated  protein  (Strap),  mRNA 


I 


V 


TC208960  2.779201 

TC174179  2.777212 

TCI  88305  2.769935 

TC178738  2.715629 

C78068  2.696100 

C87638  2.679857 

TCI  57385  2.675020 

AU041329  2.651487 

L0208A04-3  2.645424 

AU021860  2.603490 

AU022194  2.598836 

NM_007561  2.594891 

TC173759  2.577267 

AU022374  2.566131 

TC174791  2.531756 

C80763  2.498239 

TCI  73399  2.490726 

H3126E09  2.477740 

AU021733  2.463330 

TC162016  2.459140 

C78061  2.443694 

C86825  2.440315 

C78093  2.437616 

AU022156  2.435370 

TCI  62632  2.404160 

AU022430  2.385197 

C78157  2.374702 

EST03029  2.366326 

TCI  89251  2.343817 

C86821  2.323865 

TCI  73324  2.310638 

TC205126  2.306227 

AU022065  2.303341 

AU021745  2.299328 

AU023256  2.296894 

TCI 92702  2.288848 

C79008  2.271746 

AU022321  2.270753 

EST03446  2.263182 

C76575  2.262311 

AU022331  2.262090 

AU024767  2.257646 

TCI  98946  2.241669 

TC187710  2.233919 

TCI  72405  2.224768 

TCI  88356  2.222745 

TCI  86921  2.220592 

EST03496  2.214832 


Mus  musculus  nuclear  transcription  factor  RelA  (Rela)  gene,  complete  cds 

Mus  musculus  Ena-VASP-like  isoform  (Evll)mRNA,  complete  cds, alternatively  spliced 

Huamn  BMP2/4 

Mus  musculus  X  transporter  protein  3  (Xtrp3),  mRNA 
Mus  musculus  signal  sequence  receptor,  delta  (Ssr4),  mRNA 

Mus  musculus  X-linked  lymphocyte-regulated  3b  (Xlr3b),  mRNA 

Rattus  norvegicus  bHLH  transcription  factor  Misti  (Misti)  gene,  complete  cds 

Homo  sapiens  KIAA0553  protein  gene,  complete  cds;and  alphallb  protein  gene 
Psammomys  obesus  beacon  mRNA,  complete  cds 


Rattus  norvegicus  mRNA  for  dihydrolipoamide  acetyltransferase 

Mus  musculus  mRNA  for  PC3B  protein 

Mus  musculus  mRNA  for  hypothetical  protein  expressed  in  thymocytes),  partial 

Mus  musculus  U22  snoRNA  host  gene  (UHG)  gene,  complete  sequence 


TCI 74040 

2.198529 

TC203464 

2.194281 

Mus  musculus  caspase  12  (Casp12),  mRNA 

L0219C11-3 

2.191261 

TC163941 

2.189350 

TCI 65384 

2.187309 

Homo  sapiens  sciellin  (SCEL),  mRNA 

AU022363 

2.185691 

C87602 

2.178928 

AU024748 

2.175758 

AU021952 

2.163344 

C80899 

2.161401 

TC200616 

2.160273 

Mus  musculus  gene  rich  cluster,  C8  gene  (Grcc8),  mRNA 

NM_0 19305 

2.156498 

Human  FGF2 

AU024596 

2.154303 

TCI  63888 

2.146388 

Mouse  integrin  beta  4  subunit  mRNA 

AU023218 

2.144509 

EST02095 

2.142997 

AU023189 

2.140748 

NM_013414 

2.139277 

Rat  OSTEOCALCIN 

AU024605 

2.137971 

TCI  65048 

2.137789 

TCI 60424 

2.136576 

Homo  sapiens  E-1  enzyme  (MASA),  mRNA 

C78176 

2.112627 

AU043475 

2.109897 

AU024499 

2.106700 

AU024601 

2.104568 

TCI  88478 

2.103129 

AU021766 

2.096848 

AU022382 

2.096010 

TC202511 

2.088855 

Mus  musculus  Pro-rich,  PH,  SH2  domain-containing  signaling  mediator  (PSM)  mRNA 

AU022202 

2.086047 

TC166103 

2.081321 

TC202623 

2.079226 

C86859 

2.078406 

Homo  sapiens,  clone  hRPK.58_A_1,  complete  sequence 

TCI 661 62 

2.067844 

EST03382 

2.067519 

Mus  musculus  transthyretin  (Ttr),  mRNA 

TC208089 

2.065965 

TC161597 

2.065945 

TC205477 

2.064887 

Mus  musculus  potassium  voltage-gated  channel, subfamily  H(eag-related),(Kcnh1) 

TC186112 

2.063845 

TC204141 

2.062597 

Rattus  norvegicus  mevalonate  pyrophosphate  decarboxylase  mRNA,  complete  cds 

TCI  97087 

2.057604 

AU022204 

2.056734 

AU022409 

2.055760 

TCI  92594 

2.053742 

EST03436 

2.046538 

TC202525 

2.043880 

Rattus  norvegicus  hypertension-related  calcium-regulated  gene  mRNA,  complete  cds 

TC2 10779 

2.042810 

Homo  sapiens  KIAA0513  gene  product  (KIAA0513),  mRNA 

AU022341 

2.042184 

AU022070 

2.042046 

AU024727 

2.041442 

AU043452 

2.036570 

TCI  88626 

2.035291 

M.musculus  ml 

EST03723 

2.032456 

AU022460 

2.029631 

TCI  60632 

2.028687 

Mus  musculus 

Z46629.1 

2.027188 

Rat  SOX  9 

TC205634 

2.024558 

Mus  musculus 

TC201670 

2.024546 

Mus  musculus 

AU022330 

2.024488 

TC161617 

2.022511 

Mus  musculus 

AU024594 

2.021603 

Mus  musculus 

C81146 

2.016730 

TC203654 

2.008460 

TC161681 

2.007609 

C87514 

2.003433 

TCI  88237 

2.002615 

Mus  musculus 

TCI 74887 

2.000214 

Mus  musculus 

M.musculus  mRNA  for  arachidonate  epidermis-type  12(S)-lipoxygenase 


As  shown  in  Table  3,  there  are  142  genes  that  are  upregulated  in  cells  subjected  to  fluid 
flow  shear  stress  as  compared  to  the  non-stressed  cells  in  C3H  mice.  There  are  134  genes  that  are 
downregulated  in  stressed  vs.  non-stressed  cells  in  these  mice  (Table  4). 


Comparison  of  gene  expression  profiles  between  B6  and  C3H  mice. 

Since  our  main  objective  is  to  find  out  the  genes  involved  in  the  mechanical  signaling,  we 
compared  the  gene  expression  profiles  between  B6  and  C3H.  A  combined  list  of  B6S/B6C  and 
C3HS/C3HC  was  made  to  compare  B6S  with  C3HS.  A  comparison  of  stress  vs.  control 
osteoblasts  showed  that  41  genes  or  ESTs  displayed  more  than  two-fold  increase  in  expression 
(Table  5)  whereas  only  1 6  genes  were  found  to  be  downregulated  (Table  6). 

Table  5.  List  of  genes  in  which  expression  level  increased  more  than  two-fold  after 
application  of  shear  stress  in  osteoblasts  isolated  from  B6  mice  as  compared  to 
osteoblasts  isolated  from  C3H  mice.  The  ratios  for  the  microarray  were  determined  by 
using  Genespring  software. 


Fold 

Gene 

Accession  # 

t  Change 

Name 

C81354 

21.9628 

C80763 

4.20637 

AF019048 

3.61608 

Mouse  RANKL 

E01057 

3.34883 

Mouse  interleukin-1  (IL-1)  precursor 

TCI  57862 

3.08421 

Homo  sapiens  mRNA;  cDNA  DKFZp586C1620  (from  clone  DKFZp586C1620);  partial  cds 

J00370 

3.03069 

Mouse  c-fos  gene;  cellular  homolog  to  viral  oncogene 

AF075717 

2.88893 

Transforming  growth  factor  beta  1 -induced  factor  2 

AF053713 

2.83105 

Mouse  Osteoprotegerin  ligand 

LI 2447 

2.8031 

Insuline-like  growth  factor  binding  protein  5  (IGFBP5) 

TCI  8801 8 

2.75545 

TC200632 

2.75369 

AB009993 

2.74642 

X76401 

2.54899 

TC202099 

2.51276 

D31942 

2.47433 

AF1 20489 

2.45316 

X58636.1 

2.41252 

TC201479 

2.39966 

J04069 

2.39514 

V00727 

2.38357 

AF068615 

2.38111 

TC162168 

2.37976 

X56848 

2.3708 

XI 5848 

2.34551 

U78048 

2.27706 

AJ009862 

2.27282 

L28108 

2.26984 

U51001 

2.25219 

Y15800 

2.23637 

E03515 

2.21036 

SM95800 

2.20214 

U64331 

2.18271 

M29464 

2.17397 

U36384 

2.16404 

AF067191 

2.14363 

NM_0 10703 

2.09572 

X53802 

2.07883 

J05265 

2.07351 

E04743 

2.03421 

X57413 

2.01756 

SU78076 

2.00177 

Homo  sapiens  CCR4-NOT  transcription  complex,  subunit  8  (CNOT8),  mRNA 
Collagen  A1(V) 

Tumor  necrosis  factor  receptor  2 
Mouse  Oncostatin  M 

Growth  hormone  receptor/growth  hormone-binding  protein 
(LEF1)  lymphoid  enhancer  binding  factor  1 

Homo  sapiens  MEBP-1  mRNA  for  MAPK-ERK  binding  protein-1,  complete  cds 
Insulin-like  growth  factor  II  (IGF-II) 
c-fos  oncogene 

Ciliary  neutotrophic  factor  receptor  alpha  precursor 
Homo  sapiens  CGI-147  protein  (LOC51651),  mRNA 
Bone  morphogenetic  protein  4  (BMP-4) 

Retinoic  acid  receptor  gamma 

Bone  morphogenetic  protein  type  II  receptor  BRK-3 

Transforming  growth  factor-beta  1 

Mouse  PTH/PTHrP 

DLX-1  gene 

Mouse  mRNA  for  G-protein  coupled  receptor  kinase  6-D 
Mouse  interleukin  6  receptor  protein 
Mouse  Myogenin  gene 
Mouse  Osteoprotegerin  (OPG) 

Mouse  platelet-derived  growth  factor  a  chain  (PDGA) 

Mouse  twist-related  bHLH  protein  Dermo-1 
Fibroblast  growth  factor  l.a  gene,  exon  1 
Lef  1 

Mouse  mRNA  for  interleukin-6  receptor 
Interferon  gamma  receptor 
Mouse  IL-1  alpha 
Transforming  growth  factor-beta  2 
Mouse  sepiaterin  reductase  gene 


Table  6.  List  of  genes  in  which  expression  level  decreased  more  than  two-fold  after 
application  of  shear  stress  in  osteoblasts  isolated  from  B6  mice  as  compared  to 
osteoblasts  isolated  from  C3H  mice.  The  ratios  for  the  microarray  were  determined  by 
using  Genespring  software. 


Fold 


Accession  # 

Change  Gene  Name 

C86942 

4.135808 

C78087 

2.929459 

EST03087 

2.875889 

TCI  78738 

2.715629 

C78068 

2.696100 

AU022194 

2.598836 

C78157 

2.374702 

TCI  92702 

2.288848 

C79008 

2.271746 

AU022331 

2.262090 

AU024767 

2.257646 

TCI  74040 

2.198529 

TCI  63941 

2.189350 

TCI  61 597 

2.065945 

TC2 10779 

2.042810 

C81146 

2.016730 

Mus  musculus  mRNA  for  PC3B  protein 


Homo  sapiens  KIAA0513  gene  product  (KIAA0513),  mRNA 


Specific  Objective  3:  Validation  of  microarray  data  with  Real  Time  PCR. 

To  further  confirm  our  microarray  data,  we  are  now  in  the  process  of  selecting  genes  for 
doing  Real  Time  PCR.  The  selection  of  genes  was  based  on:  1)  confirmation  of  the  genes  known 
to  be  involved  in  the  mechanical  signaling;  and  2)  finding  the  genes  or  ESTs  with  unknown 
function  that  may  be  involved  in  the  mechanical  signaling.  Thus,  Real  Time  PCR  was  performed 
for  a  known  gene  (c-fos)  and  two  ESTs  (  C81354  and  C80763)  selected  from  the  microarray 
gene  list.  The  results  show  that  shear  stress  induces  upregulation  of  c-fos  and  C80763  in  B6  bone 
cells  and  not  in  C3H  bone  cells.  These  results  were  consistent  with  the  data  derived  from  the 
microarray  analysis  (Table  7).  However,  in  case  of  C81354,  microarray  results  showed  a  21 -fold 
increase  in  expression  in  B6  mice  as  compared  to  C3H  mice,  whereas  real  time  PCR  showed  no 
difference  in  the  expression  level  in  B6  and  C3H  osteoblasts  This  shows  a  discrepancy  between 
the  microarray  and  real  time  PCR  data  and  emphasizes  the  importance  of  confirming  the 
microarray  data  w  ith  other,  more  authentic  techniques  as  real  time  PCR. 

Above  results  confirm  that  c-fos  is  involved  in  the  mechanical  signaling  and  further 
suggest  that  EST  (accession  #  C80763)  might  play  an  important  role  in  mechanical  stress  induced 
cell  proliferation  and  differentiation  which  emphasizes  the  importance  of  studying  this  EST  for  its 
role  in  mechanical  signaling. 


Gene  Fold  Change 

(Stress  vs.  Control) 

_ Microarrav _ Real  Time  PCR 


EST  (accession# 

4.20 

4.18 

C80763 

EST  (accession  # 

21.96 

1.67 

C81354 

c-fos 

3.03 

2.92 

Table  7.  Gene  expression  change  after  fluid  flow  shear  stress  in  osteoblasts  isolated  from 
B6  mice  by  microarray  and  Real  Time  PCR.  The  fold  changes  were  defined  in  relative  to 
the  expression  level  of  osteoblasts  not  subjected  to  shear  stress. 

Additional  Progress  for  Technical  Objective  3: 


In  addition  to  C3H  and  B6  mice  strains,  we  evaluated  bone  formation  response  to  loading 
in  6  other  inbred  strains  of  mice  (129J,  BALB/cByJ,  NZB/BINJ,  RF/J,  AKR/J,  and  CBA/J). 
These  additional  strains  of  mice  were  selected  for  the  following  reasons:  1)  to  establish  the  mouse 
strain  that  shows  the  best  response  to  loading;  and  2)  to  select  an  optimal  mouse  pair  for  studies 
on  evaluation  of  candidate  genes  involved  in  mediating  loading  response  on  bone.  Osteoblasts 
isolated  from  these  mice  strains  were  subjected  to  fluid  flow  shear  strain  and  thymidine 
incorporation  was  measured  after  24  hours  of  stress.  As  shown  in  Figure  6,  CBA/J,  B6,  RF/J, 
Balbc  and  NZB  show  a  significant  increase  in  the  shear  stress  induced  cell  proliferation  whereas 
C3H,  129/J  and  AKR  showed  no  significant  change  in  the  shear  stress  induced  cell  proliferation 
as  compared  to  the  normal  cells  from  the  respective  strains  of  mice. 
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Figure  6:  Effect  of  fluid  flow  shear  strain  on  [3H]thymidine  incorporation  in  bone  cells  isolated  from 
different  strains  of  mice.  Cells  were  subjected  to  fluid  flow  for  30  min  and  [3H]Thymidine  incorporation 
was  measured  after  24  hours  of  incubation.  **p<0.01,  *p<0.05. 

Based  on  the  response  to  in  vitro  mechanical  loading  (shear  stress),  the  overall  ranking  of 
different  mouse  strains  is  shown  in  Table  8.  Accordingly,  C3H  and  129/J  seem  to  be  poor 
responders  whereas  B6,  RF/J,  Balbc,  and  NZB  seem  to  be  good  responders. 


Rank 

Mouse  Strain 

Thymidine  Incorporation 
(%  of  control) 

Response  to 
Shear  Stress 

1 

C3H 

108  ±22.69 

Poor 

2 

129J 

110  ±08.46 

Poor 

3 

AKR 

142  ±  14.85 

Medium 

4 

CBA 

145  ±  47.06 

Medium 

5 

B6 

183  ±46.38 

Good 

6 

RFJ 

190  ±24.70 

Good 

7. 

Balbc 

195  ±  12.93 

Good 

8 

NZB 

208  ±  70.58 

Good 

Table  8.  Ranking  of  mouse  strains  (Good  and  Poor  responder)  based  upon  the  thymidine 
incorporation  in  response  to  fluid  flow  shear  stress  in  osteoblasts  isolated  from  different 
mouse  strains. 

Reportable  Outcomes 
Publication: 

Sonia  Kapur,  David  J.  Baylink,  and  K.-H.  William  Lau:  Fluid flow  shear  stress 
stimulates  human  osteoblast  proliferation  and  differentiation  through  multiple  interacting 
and  competing  signal  transduction  pathways.  Bone  (2003);  32:  241-251. 

Abstract: 

Sonia  Kapur,  S.T  Chen,  D.  J.  Baylink,  K.-H.  W.  Lau:  Extracellular  Signal-Regulated 
Kinase  (ERK)-1  and  ERK2  Are  Both  Essential  in  the  Mediation  of  The  Flow  Shear 
Strain-Induced  Human  Osteoblast  Proliferation.  (2002);  At  24th  Annual  Meeting  of 
American  Society  for  Bone  and  Mineral  Research.  San  Antonio,  Texas.  September  20- 
24. 

Key  Findings 

1.  We  evaluated  strain-induced  tyrosine  phosphorylation  levels  of  key  signaling  proteins  in 
C3H  and  B6  mice.  We  studied  the  fluid  flow  induced  phosphorylation  of  MAPK  and 
integrin  expression  in  these  cells  and  found  that  shear  strain  induced  a  significant 
increase  in  the  integrin  plexpression  as  well  as  the  phosphorylation  levels  of  ERKl  and 
ERK2  in  B6  bone  cells.  On  the  contrary,  no  change  was  observed  in  either  integrin 
Plexpression  or  phosphorylation  levels  of  ERK1  and  ERK2  in  C3H  bone  cells  in 
response  to  fluid  flow. 

2.  We  also  found  that  the  protein  kinase  activity  of  both  ERK1  and  ERK2  are  essential  for 
the  fluid  flow  shear  stress-induced  osteoblast  proliferation. 

3.  The  most  intriguing  finding  of  this  study  is  that  overexpressing  of  an  inactive  form  of 
ERK1  or  ERK2  alone  is  sufficient  to  block  the  actions  of  endogenous  ERKs  in  response 
to  flow  shear  stress  in  bone  cell  proliferation,  further  suggesting  that  both  ERK1  and 
ERK2  are  essential  for  shear  stress-induced  proliferation  of  bone  cells. 

4.  The  microarray  analysis  revealed  that  with  over  5,000  genes  explored,  109  genes  show 
more  than  a  two-fold  increase  and  139  genes  show  more  than  two-fold  decrease  in 
expression  in  cells  subjected  to  fluid  flow  shear  stress  as  compared  to  the  non-stressed 
cells  in  B6  mice.  In  case  of  C3H  mice,  there  are  142  genes  that  show  more  than  two-fold 
increase  in  expression  and  135  genes  that  are  downregulated  more  than  two-fold  in  cells 
subjected  to  fluid  flow  shear  stress  as  compared  to  the  non-stressed  cells.  A  comparison 
of  stress  vs.  control  osteoblasts  showed  that  only  41  genes  or  ESTs  displayed  more  than 
two-fold  increase  whereas  only  16  genes  or  ESTs  showed  more  than  two-fold  decrease  in 
expression  between  B6  and  C3H  mice. 

5.  Microarray  and  real  time  PCR  results  confirmed  that  c-fos  is  involved  in  the  mechanical 
signaling  and  further  suggested  that  EST  (accession  #  C80763)  might  play  an  important 
role  in  mechanical  stress  induced  cell  proliferation.  In  future  studies,  we  shall 
characterize  C80763  and  study  its  role  in  mechanical  signaling. 

6.  In  addition  to  B6  and  C3H,  bone  formation  response  to  shear  stress  in  six  different  inbred 
strains  of  mice  showed  that  RF/J,  Balbc  and  NZB  are  responsive  to  shear  stress  whereas 
129/J  is  non-responsive  to  shear  stress.  We  are  now  in  the  process  of  determining 


whether  some  combination  of  these  six  strains  would  be  more  appropriate  to  the  C3H  and 
B6  for  the  high  and  low  response  pair. 

Conclusions 

1.  The  in  vitro  studies  show  that  shear  stress  induces  an  increase  in  the  phosphorylation 
levels  of  ERK1  and  ERK2  in  B6  bone  cells  whereas  no  change  was  seen  in  C3H  cells. 
Our  studies  further  show  that  the  protein  kinase  activity  of  both  ERK1  and  ERK2  are 
essential  for  the  fluid  flow  shear  stress-induced  osteoblast  proliferation. 

2.  The  microarray  data  reveals  the  list  of  candidate  genes  responsible  for  mediating  the  bone 
formation  response  to  mechanical  loading.  With  over  5,000  genes  explored,  41  of  these 
genes  showed  more  than  two-fold  increase  in  expression  in  B6  bone  cells  above  C3H 
cells  when  subjected  to  shear  stress.  These  genes  will  now  be  used  to  further  evaluate  the 
mechanical  stress  differential  responses  between  these  two  inbred  strains  of  mice. 

3.  Based  on  in  vitro  studies  on  bone  formation  response  to  shear  stress  in  eight  different 
strains  of  mice,  C3H  and  B6  are  appropriate  for  the  poor  and  good  response  pair  to 
mechanical  loading. 
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2.  Genetic  Analysis  of  Three  Large  Pedigrees  with  Very  High  Bone  Density 
Introduction 

Quantitative  traits,  such  as  bone  mineral  density  (BMD),  are  usually  due  to  genetic 
contributions  from  multiple  genes.  Our  long-term  goal  in  this  project  is  to  identify  those  genes 
that  contribute  to  peak  bone  density  by  linkage  analysis  studies  in  families  with  very  high  BMD 
(>2.5  SD  of  normal). 

As  reported  in  last  year’s  progress  report,  three  large  pedigrees  with  affected  members 
exhibiting  bone  density  of  more  than  2.5  SD  above  the  normal  mean  were  collected.  Blood  was 
drawn  from  the  family  members  in  Argentina,  frozen  and  sent  to  the  Musculoskeletal  Disease 
Center  (MDC)  for  processing  and  analysis.  At  the  MDC  last  year,  genomic  DNA  was  purified 
and  serum  marker  levels  of  bone  formation  and  resorption  were  measured.  The  second  year  of 
the  project  has  focused  on  high-throughput  genotyping  for  gene  discovery  through  linkage 
analysis.  Bone  density  is  an  important  determinant  of  susceptibility  to  fracture.  Thus,  the 
identification  of  the  genes  that  contribute  to  high  bone  density  in  these  families  may  represent  a 
major  advance  in  understanding  the  pathways  that  regulate  bone  formation  and  bone  healing  and 
the  pathogenesis  of  bone  diseases  such  as  osteoporosis. 

A.  Technical  Objectives: 

Following  are  the  specific  objectives  during  the  second  year  of  the  grant  period: 

1)  We  will  implement  a  whole  genome  screen  with  the  ABI  377  system  using 
multiplex  PCR  amplification  of  primer  sets  labeled  with  fluorescent  labels. 

2)  Products  will  be  analyzed  by  electrophoresis  in  the  ABI  377  and  alleles  scored 
with  appropriate  software. 

3)  We  will  perform  linkage  analysis  using  the  CRI-MAP,  GENEHUNTER,  and 
FASTLINK  4.0  software  programs. 

4)  At  chromosomal  locations  showing  evidence  of  linkage  in  any  of  the  families, 
additional  informative  markers  will  be  genotyped  to  further  narrow  down  the 
region  of  the  genetic  locus. 

We  have  accomplished  all  of  the  above  specific  objectives.  Our  progress  in  each  of  the 
Specific  Objectives  is  given  below. 


Progress  on  Technical  Objectives 

Specific  Objective  1:  A  whole  genome  screen  was  implemented  and  completed  for  Argentina 
family  E.  400  fluorescently  labeled  PCR  reactions  spaced  at  a  genetic  density  of  10 
centimorgans  were  amplified  for  all  68  individuals  of  the  E  family.  Since  a  centimorgan 
corresponds  to  roughly  one  million  basepairs,  this  corresponds  to  approximately  one 
microsatellite  marker  every  10,000,000  basepairs  to  cover  the  entire  human  genome.  Table  1  lists 
the  microsatellite  markers  used  in  this  study. 

Genomic  DNA  isolated  from  these  families  is  a  precious  resource  that  is  difficult  to 
obtain.  Thus,  our  PCR  reactions  were  optimized  to  use  the  minimal  amount  of  genomic  DNA 
while  still  amplifying  robustly  and  uniformly.  The  optimal  PCR  reaction  conditions  we 
developed  are  as  follows: 


Reaction  Mixes 


9  ul  of  True  Allele  Mix  (Applied  Biosystems) 
1  ul  of  Primer  mix  (5  uM  each  primer) 

5  ul  of  genomic  DNA  (2  ng/uB _ 

Total  volume  =  15  ul 


Thermal  Cycling 


Step  # 

Cycles 

Temperature 

Time 

1 

1 

95°C 

12:00 

2 

45 

95°C 

00:10 

55°C 

00:15 

72°C 

00:30 

3 

1 

72°C 

30:00 

4 

1 

4°C 

infinity 

Table  1.  Microsatellite  markers  used  in  this  study  and  their  genetic  location. 


Chromosome 

Marker  Name 

Genetic  Distance  (centimorgans) 

1 

D1S468 

0 

1 

D1S214 

9.4 

1 

D1S450 

15.9 

1 

D1S2667 

19.9 

1 

D1S2697 

32.9 

1 

D1S199 

40.7 

1 

D1S234 

50.2 

1 

D1S255 

60.5 

1 

D1S2797 

68 

1 

D1S2890 

79.4 

1 

D1S230 

90 

1 

D1S2841 

102 

1 

D1S207 

110 

1 

D1S2868 

124.3 

1 

D1S206 

133.6 

1 

D1S2726 

142.6 

1 

D1S252 

150.27 

1 

D1S498 

155.89 

1 

D1S484 

167.3 

1 

D1S2878 

176.3 

1 

D1S196 

181.3 

1 

D1S218 

190.2 

1 

D1S238 

201.1 

1 

D1S413 

209.6 

1 

D1S249 

219.4 

1 

D1S425 

230.8 

1 

D1S213 

241 

1 

D1S2800 

250 

1 

D1S2785 

264 

1 

D1S2842 

273.1 

1 

D1S2836 

284.4 

2 

D2S319 

7.6 

2 

D2S2211 

15.61 

2 

D2S162 

20.03 

2 

D2S168 

27.06 

2 

D2S305 

38.87 

2 

D2S165 

47.43 

2 

D2S367 

54.96 

2 

D2S2259 

64.29 

2 

D2S391 

70.31 

2 

D2S337 

80.69 

2 

D2S2368 

85.48 

2 

D2S286 

94.05 

2 

D2S2333 

103.16 

2 

D2S2216 

111.21 

2 

D2S160 

122.96 

2 

D2S347 

131.51 

2 

D2S112 

141.62 

2 

D2S151 

152.04 

2 

D2S142 

161.26 

2 

D2S2330 

169.41 

2 

D2S335 

175.91 

2 

D2S364 

186.21 

2 

D2S117 

194.45 

2 

D2S325 

204.53 

2 

D2S2382 

213.49 

2 

D2S126 

221.13 

2 

D2S396 

232.9 

2 

D2S206 

240.79 

2 

D2S338 

250.54 

2 

D2S125 

260.63 

3 

D3S1297 

8.31 

3 

D3S1304 

22.33 

3 

D3S1263 

36.1 

3 

D3S2338 

42.1 

3 

D3S1266 

52.6 

3 

D3S1277 

61.52 

3 

D3S1289 

71.41 

3 

D3S1300 

80.32 

3 

D3S1285 

91.18 

3 

D3S1566 

97.75 

3 

D3S3681 

109.22 

3 

D3S1271 

117.76 

3 

D3S1278 

129.73 

3 

D3S1267 

139.12 

3 

D3S1292 

146.6 

3 

D3S1569 

158.38 

3 

D3S1279 

169.6 

D3S1614 

177.75 

D3S1565 

186.04 

D3S1262 

201.25 

D3S1580 

207.73 

D3S1311 

224.88 

D4S412 

4.74 

D4S2935 

13.96 

D4S403 

25.9 

D4S419 

33.42 

D4S391 

43.59 

D4S405 

56.95 

D4S1592 

69.53 

D4S392 

78.97 

D4S2964 

88.35 

D4S1534 

95.09 

D4S414 

100.75 

D4S1572 

107.95 

D4S406 

117.06 

D4S402 

124.45 

D4S1575 

132.05 

D4S424 

144.56 

D4S413 

157.99 

D4S1597 

169.42 

D4S1539 

176.19 

D4S415 

181.36 

D4S1535 

195.06 

D4S426 

206.98 

D5S1981 

1.72 

D5S406 

11.85 

D5S630 

19.67 

D5S419 

39.99 

D5S426 

51.99 

D5S418 

58.55 

D5S407 

64.67 

D5S647 

74.07 

D5S424 

81.95 

D5S641 

92.38 

D5S428 

95.4 

D5S644 

104.76 

D5S433 

111.97 

D5S2027 

119.5 

D5S471 

129.83 

D5S2115 

138.64 

D5S436 

147.49 

D5S410 

156.47 

D5S422 

164.19 

D5S400 

174.81 

D5S408 

195.49 

D6S1574 

9.18 

D6S309 

14.07 

6 

D6S470 

18.22 

6 

D6S289 

29.93 

6 

D6S422 

35.66 

6 

D6S276 

44.41 

6 

D6S1610 

53.81 

6 

D6S257 

79.92 

6 

D6S460 

89.83 

6 

D6S462 

99.01 

6 

D6S434 

109.19 

6 

D6S287 

121.97 

6 

D6S262 

130 

6 

D6S292 

136.97 

6 

D6S308 

144.46 

6 

D6S441 

154.1 

6 

D6S1581 

164.78 

6 

D6S264 

179.07 

6 

D6S446 

189 

6 

D6S281 

190.14 

7 

D7S531 

5.28 

7 

D7S517 

7.44 

7 

D7S513 

17.74 

7 

D7S507 

28.74 

7 

D7S493 

34.69 

7 

D7S516 

41.69 

7 

D7S484 

53.5 

7 

D7S510 

59.93 

7 

D7S519 

69.03 

7 

D7S502 

78.65 

7 

D7S669 

90.42 

7 

D7S630 

98.44 

7 

D7S657 

104.86 

7 

D7S515 

112.32 

7 

D7S486 

124.08 

7 

D7S640 

137.83 

7 

D7S684 

147.22 

7 

D7S661 

155.1 

7 

D7S636 

162.33 

7 

D7S798 

168.98 

7 

D7S2465 

180.24 

8 

D8S264 

0.73 

8 

D8S277 

8.34 

8 

D8S550 

21.33 

8 

D8S258 

41.55 

8 

D8S1771 

50.05 

8 

D8S505 

60.87 

8 

D8S285 

71 

8 

D8S260 

79.36 

8 

D8S270 

103.69 

8 

D8S1784 

118.15 

8 

D8S514 

130 

8 

D8S284 

143.82 

8 

D8S272 

154.02 

9 

D9S288 

9.83 

9 

D9S286 

18.06 

9 

D9S285 

29.52 

9 

D9S157 

32.24 

9 

D9S171 

42.73 

9 

D9S161 

51.81 

9 

D9S1817 

59.34 

9 

D9S273 

65.79 

9 

D9S175 

70.33 

9 

D9S167 

83.41 

9 

D9S283 

94.85 

9 

D9S287 

103.42 

9 

D9S1690 

106.63 

9 

D9S1677 

117.37 

9 

D9S1776 

123.33 

9 

D9S1682 

132.09 

9 

D9S290 

140.86 

9 

D9S164 

147.91 

9 

D9S1826 

159.61 

9 

D9S158 

161.71 

10 

D10S249 

2.13 

10 

D10S591 

13.49 

10 

D10S189 

19 

10 

D10S547 

29.15 

10 

D10S1653 

40.36 

10 

D10S548 

45.85 

10 

D 1  OS 1 97 

52.1 

10 

D10S208 

60.64 

10 

D10S196 

70.23 

10 

D10S1652 

80.77 

10 

D10S537 

91.13 

10 

D10S1686 

105.04 

10 

D10S185 

116.34 

10 

D10S192 

124.27 

10 

D10S597 

128.73 

10 

D10S1693 

137.39 

10 

D10S587 

147.57 

10 

D10S217 

157.89 

10 

D10S1651 

168.77 

10 

D10S212 

170.94 

11 

D11S4046 

2.79 

11 

D11S1338 

12.92 

11 

D11S902 

21.47 

11 

DUS904 

33.57 

11 

D11S935 

45.94 

11 

D11S905 

51.95 

11 

Dll S4 191 

60.08 

11 

D11S987 

67.48 

11 

D11S1314 

73.64 

11 

D11S937 

79.98 

11 

D11S901 

85.48 

11 

D11S4175 

91.47 

11 

D11S898 

98.98 

11 

D11S908 

108.59 

11 

D11S925 

118.47 

11 

Dll S4 151 

127.33 

11 

D1 1 SI 320 

141.91 

11 

D11S968 

147.77 

12 

D12S352 

0 

12 

D12S78 

11.87 

12 

D12S99 

12.6 

12 

D12S336 

19.68 

12 

D12S364 

30.6 

12 

D12S310 

36.06 

12 

D12S1617 

44.03 

12 

D12S345 

53.09 

12 

D12S85 

61.34 

12 

D12S368 

66.03 

12 

D12S83 

75.17 

12 

D12S326 

86.4 

12 

D12S351 

95.56 

12 

D12S346 

104.65 

12 

D12S79 

125.31 

12 

D12S86 

134.54 

12 

D12S324 

147.17 

12 

D12S1659 

155.94 

12 

D12S1723 

164.63 

13 

D13S175 

6.03 

13 

D13S217 

17.21 

13 

D13S171 

25.08 

13 

D13S218 

32.9 

13 

D13S263 

38.32 

13 

D13S153 

45.55 

13 

D13S156 

55.85 

13 

D13S170 

63.9 

13 

D13S265 

68.73 

13 

D13S159 

79.49 

13 

D13S158 

84.87 

13 

D13S173 

93.52 

13 

D13S1265 

98.82 

13 

D13S285 

110.55 

14 

D14S261 

6.46 

14 

D14S283 

13.89 

14 

D14S275 

28.01 

14 

D14S70 

40.11 

14 

D14S288 

47.51 

14 

D14S276 

56.36 

14 

D14S63 

69.18 

14 

D14S258 

76.28 

14 

D14S74 

87.36 

14 

D14S68 

95.9 

14 

D14S280 

105 

14 

D14S65 

117.3 

14 

D14S985 

126.61 

14 

D14S292 

134.3 

15 

D15S1002 

14.58 

15 

D15S165 

20.24 

15 

D15S1007 

25.86 

15 

D15S1012 

35.95 

15 

D15S994 

40.25 

15 

D15S978 

45.62 

15 

D15S117 

51.21 

15 

D15S153 

62.4 

15 

D15S131 

71.28 

15 

D15S205 

78.92 

15 

D15S127 

86.81 

15 

D15S130 

100.59 

15 

D15S120 

112.58 

16 

D16S423 

10.36 

16 

D16S404 

18.07 

16 

D16S3075 

23.28 

16 

D16S3103 

32.07 

16 

D16S3046 

40.65 

16 

D16S3068 

48.53 

16 

D16S3136 

62.11 

16 

D16S415 

67.62 

16 

D16S503 

83.55 

16 

D16S515 

92.1 

16 

D16S516 

100.39 

16 

D16S3091 

111.12 

16 

D16S520 

125.82 

17 

D17S849 

0.63 

17 

D17S831 

6.6 

17 

D17S938 

14.69 

17 

D17S1852 

22.24 

17 

D17S799 

31.96 

17 

D17S921 

36.14 

17 

D17S1857 

43.01 

17 

D17S798 

53.41 

17 

D17S1868 

64.16 

17 

D17S787 

74.99 

17 

D17S944 

82.56 

17 

D17S949 

93.27 

17 

D17S785 

103.53 

17 

D17S784 

116.86 

17 

D17S928 

126.46 

18 

D18S59 

0 

18 

D18S63 

8.3 

18 

D18S452 

18.7 

18 

D18S464 

31.17 

18 

D18S53 

41.24 

18 

D18S478 

52.86 

18 

D18S1102 

62.84 

18 

D18S474 

71.32 

18 

D18S64 

84.8 

18 

D18S68 

96.48 

18 

D18S61 

105.03 

18 

D18S1161 

114.26 

18 

D18S70 

126 

19 

D19S209 

10.97 

19 

D19S216 

20.01 

19 

D19S884 

26.37 

19 

D19S221 

36.22 

19 

D19S226 

42.28 

19 

D19S414 

54.01 

19 

D19S220 

62.03 

19 

D19S420 

66.3 

19 

D19S902 

72.72 

19 

D19S571 

84.08 

19 

D19S418 

92.56 

19 

D19S210 

100.01 

20 

D20S117 

2.83 

20 

D20S889 

11.2 

20 

D20S115 

21.15 

20 

D20S186 

32.3 

20 

D20S112 

39.25 

20 

D20S195 

50.81 

20 

D20S107 

55.74 

20 

D20S119 

61.77 

20 

D20S178 

66.16 

20 

D20S196 

75.01 

20 

D20S100 

84.78 

20 

D20S171 

95.7 

20 

D20S173 

98.09 

21 

D21S1256 

9.72 

21 

D21S1914 

19.39 

21 

D21S263 

27.4 

21 

D21S1252 

35.45 

21 

D21S266 

45.87 

22 

D22S420 

4.06 

22 

D22S539 

14.44 

22 

D22S315 

21.47 

22 

D22S280 

31.3 

22 

D22S283 

38.62 

22 

D22S423 

46.42 

22 

D22S274 

51.54 

Following  amplification,  the  400  markers  were  pooled  into  28  multiplex  panels  and 
electrophoresed  with  a  high-density  LIZ  labeled  size  standard.  Over  26,000  genotypes  were 
generated  for  this  family  by  this  method.  Figure  1  shows  a  multiplex  pool  of  panel  28. 


Figure  1.  Multiplex  pool  of  18  fluorescently  labeled  PCR  products  amplified  from  an  individual 
of  family  E.  All  18  reactions  were  run  in  a  single  capillary  on  the  ABI  3100  DNA  analyzer. 
PCR  reactions  are  labeled  with  FAM  (blue),  VIC  (green),  and  NED  (yellow,  shown  as  black  in 
the  plot  above  for  clarification).  The  labeled  genotyping  PCR  products  are  co-electrophoresed 
with  a  high  density  size  standard  labeled  with  LIZ  (orange).  The  x-axis  is  size  in  base  pairs  of  the 
respective  markers.  The  y-axis  is  fluorescent  intensity. 

Specific  Objective  2:  Following  electrophoresis  on  the  ABI  3 1 00  DNA  Analyzer,  Genotyper 
software  macros  were  used  to  semi-automatically  score  the  allele  calls  for  all  the  28  multiplex 
pools  of  every  E  family  individual.  After  initial  scoring  by  these  macros,  allele  calls  were 
visually  checked  and  edited  if  necessary.  A  table  of  the  calls  was  generated  and  the  allele  sizes 
converted  into  respective  allele  bins.  Figure  2  and  Table  1  demonstrates  the  method  used  to 
generate  all  26,000  allele  calls  for  family  E. 


Figure  2.  Genotyper  plots  of  allele  calls  for  individual  9  of  E  family.  The  data  shows  multiplex 
pool  28.  The  upper  blue  panel  is  6  PCR  reactions  labeled  with  FAM.  The  center  plots  are  7  PCR 
reactions  labeled  with  VIC  (green)  and  5  PCR  reactions  labeled  with  NED  (yellow  dye  but  black 
in  plot  for  clarity).  The  lowest  orange  plot  is  the  LIZ  labeled  molecular  size  standard.  The  macro 
has  labeled  the  allele  peaks  with  the  size  in  nucleotide  basepairs.  Again  the  x-axis  is  size  in 
nucleotides  and  the  y-axis  is  fluorescent  intensity. 


Table  1.  Genotyper  macro  output  table  of  E  family  individual  9  for  chromosome  1.  The 
generated  table  lists  the  person  ID,  the  capillary  number  the  sample  was  run  on  (AbilD),  the  dye 
label  color,  the  marker  name,  the  multiplex  panel  pool,  the  chromosome  the  marker  is  located  on, 
the  allele  sizes  of  the  marker,  the  calculated  allele  bin  numbers  and  whether  the  marker  is 
heterozygous  or  homozygous  for  that  individual. 
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Specific  Objective  3:  For  generation  of  the  Argentinian  pedigrees,  Progeny  software  was  used. 
Progeny  is  a  computer  program  with  an  integrated,  customizable  relational  database  for  data 
management  and  pedigree  analysis.  With  Progeny  one  can  draw  pedigree  charts,  query,  sort,  and 
import  and  export  genetic  data  using  it’s  integrated  spreadsheet  module.  Figure  3  shows  pedigree 
charts  of  the  E  family  generated  using  Progeny. 

Prior  to  linkage  analysis  and  following  generation  of  the  allele  calls  and  binning  of  the 
alleles,  all  the  data  was  tested  for  potential  genotyping  errors  and  miscalls  (i.e.  non-hereditary  due 
to  Mendel  errors)  using  Pedcheck  software  (1).  Most  errors  found  were  due  to  incorrect  binning 
of  alleles.  All  miscalls  and  binning  errors  were  corrected  and  the  data  tables  exported  as  tab 
delimited  files  for  linkage  analysis. 


aoohoo: 


Figure  3a.  Pedigree  chart  of  the  Argentina  E  family.  The  chart  was  generated  using  Progeny  Software. 
Six  generations  of  the  E  family  are  shown. 


Figure  3b.  Pedigree  chart  of  smaller  region  of  E  family  containing  many  of  the  individuals  with 
high  bone  density.  Circles  indicate  females  and  squares  males.  Individuals  with  a  slash  through 
them  are  deceased.  The  number  below  the  circle/square  indicates  the  coded  individual  identifier 
and  beneath  this  the  individual’s  age.  The  pedigree  is  color  coded  for  sample  tracking.  Blue 
indicates  that  a  DEXA  BMD  score  was  done  on  the  individual,  red  indicates  that  blood  was 
drawn  and  received  by  the  MDC,  yellow  indicates  that  the  individuals  completed  the  clinical 
questionnaire,  and  center  black  circles  indicate  that  the  individual  has  very  high  bone  density  for 
both  spine  and  hip  (>2.5  SD  above  mean). 

GENEHUNTER  was  used  for  linkage  analysis  of  the  genotyping  data  (2,  3).  Several 
regions  of  the  genome  were  found  that  give  highly  significant  LOD  scores.  A  LOD  score  is  the 
logarithm  of  odds  score,  ie.  the  logarithm  of  the  likelihood  that  two  loci  are  linked/likelihood  that 
loci  are  unlinked.  Thus  a  LOD  score  of  2  indicates  that  there  are  100:1  odds  that  a  respective 
genetic  region  shows  linkage  to  that  trait.  Figure  4  shows  LOD  plots  for  chromosomes  11  and 
1 8  for  the  affected/unaffected  GENEHUNTER  linkage  analysis  of  the  E  family.  Regions  of  the 
genome  that  give  highly  significant  LOD  scores  include  parts  of  chromosomes  11  and  18. 
Additional  evidence  of  linkage  to  these  regions  was  confirmed  by  analyzing  the  data  using 
SOLAR  linkage  analysis  software.  SOLAR  stands  for  Sequential  Oligogenic  Linkage  Analysis 
Routines,  and  was  developed  by  Blangero  and  Almasy  specifically  for  linkage  analysis  of 
quantitative  traits  in  pedigrees  of  large  size  and  complexity  (4,  5).  The  CRI-MAP  and 
FASTLINK  4.0  software  programs  mentioned  in  earlier  specific  objectives  were  developed 
nearly  ten  years  ago  and  newer  linkage  analysis  packages  have  largely  superceded  them. 


Figure  4.  LOD  Score  plots  of  GENEHUNTER  linkage  analysis  for  the  E  family  for 
chromosomes  1 1  and  1 8. 
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Chromosome  1 1  GENEHUNTER  LOD  Score  plot.  LOD  scores  of  2  are  found  at  approximately 
70  centimorgans  form  1  lptel. 
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Chromosome  18  GENEHUNTER  LOD  Score  plot.  LOD  scores  of  over  2  are  found  at 
approximately  10  centimorgans  form  18ptel. 


Specific  Objective  4:  One  of  the  highly  significant  regions  in  the  genome  that  gives  LOD 
scores  over  2  is  a  region  of  chromosome  11.  In  other  studies,  the  LRP5  gene  in  this  region  of 
chromosome  1 1  was  found  to  have  missense  mutations  that  contribute  to  high  bone  density  (6,  7, 
8).  Thus  in  all  the  Argentina  families  we  are  currently  in  the  process  of  screening  for  mutations 
in  all  exons  and  splice  junctions  of  the  LRP5  gene.  If  no  mutations  are  found,  denser  targeted 
genotyping  of  chromosomes  1 1  and  18  will  be  done  in  the  other  Argentinian  families.  Also,  we 
are  currently  examining  the  Celera  and  public  genome  databases  to  determine  functional  and 
positional  candidate  genes  in  each  of  the  regions  with  significant  LOD  scores.  This  will  be 
followed  up  by  mutation  screening  of  strong  candidate  genes  and  by  denser  more  targeted 
genotyping  of  those  regions  to  reduce  the  genetic  distance  of  regions  showing  statistically 
significant  genetic  linkage.  This  will  also  reduce  the  number  of  candidate  genes  that  will  be 
needed  to  be  screened  for  mutations  in  order  to  discover  the  gene(s)  and  mutation(s)  contributing 
to  high  bone  density. 


Additional  Progress 

Evaluation  of  Biochemical  Markers 

To  assess  if  the  increased  bone  density  in  high  bone  density  family  members  is  due  to  increased 
bone  formation  and/or  due  to  decreased  bone  resorption,  we  measured  two  bone  formation 
markers  (serum  bone  specific  alkaline  phosphatase  and  N-terminal  Type-I  procollagen  peptide) 
and  one  bone  resorption  marker  (type  I  collagen  cross-linked  N-telopeptide,  NTx).  All  serum 
samples  were  analyzed  in  one  batch  in  order  to  minimize  the  inter-assay  variation.  Before 
measurements  were  initiated,  all  assays  were  validated  for  sensitivity,  intra-  &  inter-assay 
variation,  interference  due  to  hemolyzed  or  ecteric  serum,  spike  recoveries,  linearity  and 
specificity.  A  brief  description  of  all  assays  is  provided  below. 

Alkaline  phosphatase 

A  kinetic  colorimetric  assay  was  used  to  measure  bone  specific  alkaline  phosphatase  (sALP). 
The  bone  specific  alkaline  phosphatase  is  determined  by  the  heat  inactivation  of  one  aliquot  of 
serum  sample,  standards  and  controls  at  54°C  for  12  minutes.  The  second  aliquot  is  assayed 
without  heating  and  represents  the  total  activity.  Heat  inactivated  normal  human  serum  spiked 
with  known  amounts  of  highly  purified  sALP  and  liver  alkaline  phosphatase  was  used  as 
standards  and  controls.  Sensitivity  of  total  alkaline  phosphatase  assay  was  0.6  U/L  and  average 
intra-  and  inter-assay  variations  was  CV<1 1%. 

N-Terminal  procollagen  peptide  (P1NP) 

Serum  levels  of  N-terminal  type-I  procollagen  peptide  were  was  measured  by  a 
radioimmunoassay  (RIA)  (manufactured  by  Orion  Diagnostica,  Finland  and  distributed  by 
DiaSorin,  Inc.,  Stillwater,  MN)  that  utilizes  a  polyclonal  antibody  directed  against  the  alpha  1 
chain  of  N-terminal  procollagen  peptide  of  type-I  collagen.  Levels  of  P1NP  reflect  the  formation 
of  type  I  collagen  and  as  such  are  a  marker  of  bone  formation.  Results  are  expressed  in 
micrograms  per  liter  (pg/liter).  The  sensitivity  is  2  pg/liter  and  average  intra-  and  inter-assay  CVs 
of  P1NP  radioimmunoassay  was  less  than  7%. 

N-Telopeptide  of  type-I  collagen  (NTx) 

To  assess  bone  resorption,  we  measured  serum  levels  of  type  I  collagen  cross-linked  N- 
telopeptides  using  an  enzyme  linked  immunoassay  (ELISA)  (Osteomark,  Ostex  International. 

Inc.,  Seattle,  WA,  USA)  that  utilizes  a  monoclonal  antibody  directed  against  the  cross-linking 


domain  of  N-terminal  peptide  of  type-I  collagen  in  urine.  Sensitivity  of  the  ELISA  was  5  nM 
BCE  (bone  collagen  equivalent)  and  the  linear  range  of  the  assay  was  5-40  nM  BCE.  The 
average  intra-  and  inter-assay  variation  for  the  controls  was  CV<9%. 


Summary  of  Biochemical  Markers  Results 

Bone  biochemical  marker  levels  in  Family-E  are  shown  in  Table  -  1.  Comparison  of 
bone  biochemical  markers  between  affected  individuals  in  Family-E  shows  no 
differences  in  serum  NTx  and  pro-collagen  peptide  (P1NP)  levels.  However,  the  skeletal 
alkaline  phosphatase  levels  show  significant  differences  when  males  and  females  were 
combined  and  compared  with  un-affected  family  members  (Table  -2).  Since  there  were 
only  two  affected  males,  we  did  not  perform  any  comparison  between  affected  males  and 
un-affected  family  members. 

There  are  two  types  of  bone  diseases  that  can  lead  to  high  bone  density:  osteopetrosis 
and  osteosclerosis.  Osteopetrosis  is  distinguished  by  a  low  bone  resorption  rate,  whereas 
osteosclerosis  is  distinguished  by  a  high  bone  formation  rate. 

Our  data  show  that  for  the  bone  resorption  parameter,  serum  NTx  was  not  decreased  in 
the  affected  compared  with  the  unaffected  individuals  pointing  toward  osteosclerosis  and 
high  bone  formation  as  the  cause  of  the  high  bone  density.  Moreover,  the  alkaline 
phosphatase  was  elevated  in  the  affected  individuals  with  the  high  bone  density  compared 
to  the  controls  in  the  pedigree  with  a  normal  bone  density.  In  aggregate,  these  data 
strongly  suggest  that  the  functional  abnormality  leading  to  the  high  bone  density  in  this 
pedigree  is  a  high  bone  formation  rate. 


Table  - 1.  Biochemical  markers  levels  in  Family-E 


ID 

Sex 

Age 

sNTX  (nM  BCE) 

PINP  (pg/L) 

SAPL  (U/L) 

1 

F 

56 

11 

83.9 

3.7 

9 

F 

53 

11 

72.9 

4.4 

14 

F 

58 

8.3 

93.6 

4 

16 

F 

48 

6.8 

60.6 

9 

17 

M 

58 

5.4 

37.6 

16 

19 

M 

57 

5.6 

41.9 

12 

20 

F 

57 

8.8 

55.5 

6.5 

21 

M 

59 

13 

98.2 

10 

22 

F 

52 

9.5 

91.4 

18 

23 

M 

56 

10 

47.8 

10 

24 

F 

45 

7.5 

29.8 

9.3 

27 

F 

59 

19 

78.2 

9.2 

28 

M 

62 

11 

45.8 

15 

30 

F 

59 

10 

69.6 

9.5 

31 

M 

50 

9.6 

41.9 

19 

32 

F 

51 

11 

119.7 

13 

35 

F 

68 

11 

68.2 

35 

36 

F 

46 

8.8 

104.9 

13 

41 

F 

57 

7.4 

31.3 

3.8 

42 

M 

69 

17 

84.4 

7.7 

43 

F 

61 

9.8 

47.4 

7 

45 

F 

25 

7.5 

52.5 

8.1 

48 

F 

29 

6.5 

28.7 

4.7 

49 

F 

26 

13 

41.9 

7.8 

53 

F 

38 

10 

65.9 

9.4 

54 

M 

25 

6.2 

70.7 

8.3 

55 

F 

35 

5.4 

30.4 

5.9 

56 

F 

33 

14 

98.4 

7.1 

57 

M 

20 

14 

66.4 

5.9 

60 

M 

34 

26 

63.4 

16 

62 

M 

25 

7 

52.8 

4.9 

63 

M 

43 

19 

39.0 

3.5 

64 

F 

32 

4.9 

33.9 

6.6 

65 

F 

28 

14 

49.2 

12 

66 

F 

22 

19 

58.7 

9.7 

67 

M 

16 

65.0 

11 

69 

F 

25 

11 

56.0 

7.2 

70 

M 

37 

15 

100.9 

5.3 

71 

F 

31 

6.8 

33.1 

3.7 

73 

M 

22 

25 

119.2 

17 

74 

F 

30 

29 

67.2 

9.1 

76 

F 

20 

18 

55.5 

17 

77 

F 

23 

24 

49.9 

14 

81 

M 

37 

15 

54.1 

9.7 

82 

M 

29 

25 

70.2 

11 

83 

F 

33 

17 

55.7 

6 

84 

F 

36 

10 

21.2 

3.4 

85 

M 

18 

29* 

163.9* 

37* 

86 

F 

26 

49* 

234.8* 

50* 

87 

F 

24 

74* 

631.9* 

51* 

90 

F 

42 

19 

40.3 

9 

91 

F 

40 

18 

19.4 

5.3 

95 

F 

37 

36* 

36.8 

7.4 

96 

M 

39 

11 

66.9 

13 

97 

F 

36 

7.3 

75.7 

5.4 

131 

F 

27 

23 

47.4 

3.5 

140 

M 

42 

27 

73.5 

8.8 

145 

M 

20 

88* 

285.8* 

51* 

164 

F 

34 

11 

32.2 

3.8 

168 

M 

39 

23 

72.8 

15 

172 

M 

31 

18 

34.3 

3.7 

173 

M 

27 

10 

60.1 

7.9 

WE  173 

M 

54 

13 

29.0 

9.2 

WE  174 

F 

50 

8.1 

61.7 

19 

*Due  to  very  high  values  (higher  than  clinical  reportable  range  for  these  analytes)  in 
these  patients,  we  excluded  these  samples  from  our  analysis. 


Table  -  2.  Comparison  of  skeletal  alkaline  phosphatase  between  affected  individuals  and  normal 
family  members 


Study  Groups 

n 

Age 

(Mean  ±  SD) 

sALP 

(Mean  ±  SD) 

Difference 

p-Value 

(Mann-Whitney 

TTest) 

*Affected  Male+Female 

7 

34.9+14.4 

12.6+4.4 

46% 

0.0299 

Normal  Male+Female  Family 
Members 

50 

41.1+12.8 

8.6  ±4.3 

Affected  Female 

5 

32.8+11.4 

11.1  ±4.4 

38% 

0.0914 

Normal  Female  Family 
Members 

31 

41.3+12.5 

8.0  ±4.3 

♦Marker  levels  are  affected  by  both  age  and  sex  and  hence  these  factors  should  be  taken  into  account  before  any 
comparison.  However,  because  of  only  few  affected  members  and  limited  number  of  age  and  sex  matched  family 
members,  we  have  performed  analysis  of  males  and  females  combined  together  as  well  as  separately. 


Key  Findings 

1 .  We  completed  a  whole  genome  screen  at  a  density  of  1 0  centimorgans.  This  screen 
involved  generating  over  26,000  genotypes  for  all  individuals  in  the  Argentina  E  family. 

2.  PCR  methods  and  fluorescent  multiplexing  pools  were  developed  and  optimized  to 
minimize  the  amount  of  genomic  DNA  needed  to  generate  a  whole  genome  screen.  The 
amount  of  genomic  DNA  we  now  need  for  a  10-cM  genome  screen  is  approximately  5 
micrograms  vs.  25  micrograms  in  standard  protocols.  This  conserves  the  rare  DNA 
samples  that  took  several  years  to  collect  and  isolate. 

3.  Genotyping  macros  were  developed  to  semi-automatically  call  and  generate  data  tables 
for  all  26,000  genotypes. 

4.  PedCheck  software  was  installed  and  utilized  to  insure  that  the  data  generated  was  high 
quality  with  no  Mendel  errors  and  suitable  for  linkage  analysis. 

5.  GENEHUNTER  and  SOLAR  linkage  analysis  was  implemented  and  completed  with  the 
E  family  data. 


6.  Regions  of  highly  statistically  significant  (>2  LOD  scores)  were  found  on  chromosomes 
11  and  18. 

7.  Public  and  private  databases  of  the  genomic  regions  of  significant  linkage  are  being 
searched  for  positional  and  functional  candidate  genes. 

8.  Additional  mutational  screening  of  strong  candidate  genes  is  being  undertaken  as  well  as 
denser  more  focused  genotyping  in  those  statistically  significant  genomic  regions. 

9.  All  serum  biochemical  marker  assays  have  been  completed.  From  the  serum  data,  we 
make  the  following  findings: 

a)  The  higher  sALP  values  in  affected  member  are  consistent  with  the  prediction 
that  increased  bone  formation  may  contribute  to  higher  bone  mineral  density  in 
these  patients. 

b)  Since  NTx  values  are  not  decreased  in  patients  with  high  bone  mineral  density, 
the  decreased  bone  turnover  is  less  likely  an  explanation  for  the  higher  bone 
mineral  density. 


Reportable  Outcomes 
None  yet. 


Conclusions 

There  is  strong  evidence  that  we  have  located  a  site  on  chromosome  11  that  is 
responsible  for  the  high  bone  density  phenotype  in  the  Argentina  E  family.  This  is  under 
further  investigation  through  additional  fine  genetic  mapping  and  by  mutational  screening 
of  candidate  genes.  We  have  found  additional  evidence  that  chromosome  1 8  may  harbor 
a  modifying  locus  in  this  family.  Additional  experiments  will  be  undertaken  to  confirm 
these  results  and  discover  the  gene(s)  and  mutation(s)  involved  in  the  high  bone  density 
phenotype. 
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